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ABSTRACT 

We present deep two-dimensional spectra of 22 candidate and confirmed Ly- 
man break galaxies (LBGs) at redsliifts 2 < z < 4 in the Hubble Deep Field 
(HDF) obtained at the Keck II telescope. The targets were preferentially selected 
with spatial extent and/or multiple knot morphologies, and we used slitmasks 
and individual slits tilted to optimize measurement of any spatially resolved kine- 
matics. Our sample is more than one magnitude fainter and is at higher redshift 
than the kinematic LEG targets previously studied by others. The median tar- 
get magnitude was /814 = 25.3, and total exposure times ranged from 10 to 50 
ks. We measure redshifts, some new, ranging from z = 0.2072 to z = 4.056, 
including two interlopers at 2; < 1, and resulting in a sample of 14 LBGs with a 
median redshift z = 2.424. The morphologies and kinematics of the close pairs 
and multiple knot sources in our sample are generally inconsistent with galaxy 
formation scenarios postulating that LBGs occur only at the bottom of the po- 
tential wells of massive host halos; rather, they support "collisional starburst" 
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models with significant major merger rates and a broad halo occupation distri- 
bution. For 13 LBGs with possible kinematic signatures, we estimate a simple 
dynamical mass, subject to numerous caveats and uncertainties, of the galaxies 
and/or their host dark matter halos. Dynamical mass estimates of individual 
galaxies range from 4 x lO^/i~^M0 to 1.1 x lO^^/i~^M0 and mass estimates of 
halos, based on close LBG pairs, range from < 10^°/i^^ to ~ lO^^/i^^M© with a 
median value 1 x lO^^M©. Comparison with a recent numerical galaxy formation 
model implies that indeed the pairwise velocities might not refiect true dynami- 
cal masses. We compare our dynamical mass estimates directly to stellar masses 
estimated for the same galaxies from SEDs, and find no evidence for a strong 
correlation. The diversity of morphologies and dynamics implies that LBGs rep- 
resent a broad range of galaxy or proto-galaxy types in a variety of evolutionary 
or merger stages rather than a uniform class with a narrow range of mass. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: 
high-redshift — galaxies: formation 



Introduction 



Lyman break galaxies (LBGs) are currently among our best windows on the universe 
of galaxies at redshifts z > 2. LBGs are especially useful for investigating galaxy forma- 
tion and evolution because they are relatively easy to find in large numbers (approaching 
2 arcmin^^ down to R < 25.5; Steidel et al 2003; Giavalisco et al 2002) and because they 
are relatively bright in the optical, thus permitting optical spectroscopic foUowup. Spec- 
troscopic redshifts of several thousand LBGs have been measured to d ate, with photomet 



ric redshifts of tens of thousands also avail able (jPaltani et al 1 120071 : ISteidel et al. 1 12004 



Ouchi et al. 1 l200ll:lRavindranath et allbooeh. and luminosity function s beyond z > A have 



been measured ([Kashikawa et al I l2006l : iBouwens fc lUingworth I l2006l : llwatat et al 1 120071 : 
Sawicki fc Thompson 1 120061 ). LBGs account for almost all of the star formation at 2; ~ 3 
observable in optical windows, and roughly half of the total star forma tion at those red- 



shifts once dust-obscured systems such as sub- mm galaxies are included (ISmail et al 1 12002 
Chapman et al II2005I : iGiavalisco II2002I ). 



Although a complete picture of LBGs and their relation to galaxies today remains 
elusive, many pieces of the picture are now available. Spectra of LBGs clearly show evidence 
of strong bursts of star formation, ranging from 1 — 100 yr~^, confirmed with deep radio 
and x-ray imaging with the VLA and Chandra by iReddy fc Steidel I (120041 ). who estimate 
a mean star formation rate (SFR) ~ 50 Mq yr"^ for UV-selected galaxies at 1.5 < z < 3. 
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Deep images with HST reveal small sizes ri/2 ~ Ah^^ kpc (with h = Ho/72 km s~^Mpc~^), 
high luminosities L ~ IO^^Lq, and diverse morphologies, including multiple small knots, 
diffuse halos, and asymmetrical linear features. Deep Chand ra images detect only ~ 3% of 



LBG s, implying that the AGN fraction is at most that much (ILaird et al. II2005I : iLaird et al 
2006h . 



The spectra of LBGs also show evidence of strong outflows, e.g. Lya emission lines 
red-shifted and inte rstellar absorption lines blueshifted with respect to systemic redshifts 
( Shaplev et al 2003h and a paucity of QS O H I absorption lines near LEG sightlines ( Adelberger 
20031 ). although iDesjacques et al I (120061 ) find that only smaller outflow bubbles are needed 
to explain the spectra, i.e. 0.5 Mpc comoving, implying that metals observed in the Lya 
forest are unlikely to come from LBGs. 

Measuring and understanding the dynamical state and masses of LBGs is important for 
placing them in the context of galaxy formation. A galaxy's mass affects its rate of accretion 
of new material, its ability to retain gas against the expulsive force of supernova winds, and 
possibly its eventual morphological type. The need to constrain LBG dynamics motivates 
the present work. 

In the local universe, total galaxy mass is usually m easured kinematically by deriving 
rotation curves or emission line widths for disk galaxies (ISofue fc Rubin Il200ll ). measuring 
internal velocity dispersions for elliptical and other spheroidal systems, and m easuring veloc- 



ity di stributions of satellite galaxies or other test particles for all galaxy types (jZaritsky et al 



19931 ). and then translating the observed velocity fleld into a dynamical or virial mass. The 
masses of larger systems such as galaxy groups and clusters, but not of their constituent 
galaxies, are likewise constrained by direct measurement of their velocity dispersions, as well 
as by X-ray emission and gravitational lensing. 

Three kinds of direct evidence to date have been used to study the dynamics of LBGs 
and LBG host halos: emission line widths and spatially-resolved kinematics (for total mass 
and dynamical state), spectral energy distributions (for stellar mass only), and clustering 
analysis. 

LBG mass measurements from redshifted Ha emissi on line widths and spatially-resolved 



kinematics generally fall in the range 0.5 — 25 x lO^^Mrri (IPettini et al 



2001 



Erb et al 



2006 



200 



1 



2006 ; iForster Schreiber et al II2006I : iBouche et al 1120071 : iNesvadba et al 
20091). Some signs o f rotation curves are seen ai z ^ 2 (IForster Schreiber et al II2006I . 



Forster Schreiber et al 



2009 



Lehnert et al I l2009l ) . and there is also evidence for superwinds, espec ially Lya eraission 



lines blue-shifted with respect to stella r and interstellar absorpt i on line s (IPettini et al 112001 



Erb et al II2004I : INesvadba et al 1120071 ). IForster Schreiber et al I (120091 ) found that one-third 
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of 62 galaxies at 1.3 < z < 2.6 they studied with the SINFONI integral field unit (IFU) on 
the Very Large Telescope (VLT) showed rotation-dominated kinematics, one-third are inter- 
acting or merging systems, and one-third are dominated by random motions; derived Ha dy- 
namical masses for the whole sample are 3x 10^— 3 x IO^^Mq, with median M = 2.7x lO^^Mrr 



0- 



Estimates of the stellar mass of LBGs have benefited from deep imaging with Spitzer 
Space Telescope^ s Infrared Array Camera (IRAC), including especially the GOODS and Ex- 
tended Groth Strip (EGS) fields. Stellar masses of LBGs at z > 2 generally range 
from 10^~^^Mf;) for luminosities L > L^, with typical median values M^, ~ 3 x lO^^Mfr 



{ Rigopoulou et al II2006I : iLabbe et al II2005I : IShapley et al II2005I : IPapovich et a. 



20051 ) . Lya emitting galaxies (LAE's) show similar stellar masses (ILai et al 



20011 : iHuang et al 
2007h . 



Samples of LBGs in the redshift range 2 < z < 6 are now large enough to constrain 
star formation history in LBGs: the median LBG mass is reported to be lower by around a 



factor of 10 at 2: ~ 6 th an at 2 ~ 3 (jVerma et al 1 120071 : iLabbe et al I l2006l : lYan et al 1 12006 



Eyles. L. P. et al 1120071 ). implying significant buildup of stellar mass through star formation 
and merging over that time, and mirroring the observed increase by roughly a factor of 10 
in the comoving stellar m ass density p^, from z = 3 to z = measured in the distant red 



galaxy (DRG) population (IRudnick et al II2OO6I : iFontana et al II2OO6I ) 



Finally, the observed clustering of LBGs has been used to constrain the mass of dark 
matter halos in which the galaxies presumably reside. Two-point correlation analysis reveals 
that LBGs cluster with each other with ty pical correlation length rp ^ 3 — 4if|7^ Mpc and 
implied halo masses Mha i p ~ 10^^~^^M(7^ (jGiavalisco fc Dickinson 2001 : Kashikawa et al 



20061 : lAdelberger I l2005al : lOuchi et al. II2OOII ). while still leaving the masses of individual 
LBGs within those halos uncertain. Clustering of LBGs with damped Lya QSO absorption 
line clouds ( DLAs) has also been detected, with implied DLA halo masses on the or der of 
2 X IO^^Mq jBouche fc Lowenthal 120041 : [souche et al II2OO5I : ICooke et al hoO^ . 



_LBGs are also observed to cluster on small scales: ICoUey et al I (119961 ) and lCoUey et al 



(119971 ) found that the two-point correlation function of faint objects in the HDF, consisting 
largely of LBGs at high redshift, peaks between 0'.'25 and 0'.'4 (~ 1 kpc at z > 1) with 
amplitude greater than 2, i.e. LBGs appear as multiple clumps of emission rather than 
monolithic sources. They interpret that multiplicity as evidence of starbursting regions 
within otherwise normal gas-rich galaxies. Many of the targets discussed in the current 
paper fall into that category of close pairs or multiple-knot sources. 

Galaxy formation theorists have used observations of LBGs extensively to check and 
revise their predictions and detailed descriptions of mass assembly, star formation, feed- 
back, morphology, and clustering over cosmic time. Hydrodynamic A^-body simulations 
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(|Ceverino et al II2OO9I : iFinlator et al II2OO6I : IWeinberg et al II2OO2I : IMo et al Ill999l : iNight et al 
20061 ) and semi-analytic models (SAMs; Cole et al 2000; Somerville et al 2001; Bower et al 
2006) of cold dark matter (CDM) structure and galaxy formation have generally succeeded in 
describing the observed photometric, stellar mass, size, star formation, luminosity function, 
and clustering properties of LBGs, especially by including varying amounts of dust extinc- 
tion. In particular, the maximum stellar masses of LBGs predicted in smoot hed particle 
hydro dynamics (SPH) models are M^^max = 10^^ at z = 6 and 10^^''^ at z = 3 (INight et al 



20061 ): sizes are predicted to be ri/2 ~ 0.6 — 2h^^ kpc and halo and internal velocity disper- 
sions to be ahaio ~ 180 — 290 km s~^ and amt ~ 70 — 120 km s~\ resp ectively; star formation 
rates are expected to be SFR ~ 15 — lOO Mtn yr~^ (Mo et al Ill999l) : and SFR is predicted 



to correlate only weakly with stellar mass (IWeinberg et al II2002I ) 



The available information therefore suggests that LBGs at z ~ 3 range in stellar mass 
from ~ 1 - 10 X lO^^M© with total mass from dynamical constraints M < lO^^M©, and 
that they reside in dark matter halos with mass M ~ IO^^Mq, with significant clustering. 
Their star formation rates range from 10 — 100 M© yr^^, with gas content and therefore 
star formation duty cycles and gas depletion timescales poorly constrained. Their comoving 
number density correponds roughly to that of luminous galaxies today, and their optical 
sizes and morphologies are diverse and compact. The general picture then is that LBGs are 
a heterogeneous set of actively star-forming galaxies or sub-galactic clumps, building blocks 
that will likely coalesce with other LBGs and/or non-LBG clumps to form more massive 
galaxies, groups, and clusters by z = 0. 

To address further the open question of the dynamics of Lyman break galaxies and their 
host halos and help place them with more detail in the larger context of galaxy formation 
and evolution, we have obtained spatially resolved optical (rest-UV) spectra with the 10-m 
Keck telescope of a sample of LBGs and searched for kinematic information that we can 
use to constrain individual LBG and LBG host halo dynamics directly. Our sample is more 
than one magnitude fainter and is at higher redshift than the kinematic targets previously 
studied by others. 



NIR IFU observatioris, esp ecially with adaptive optics (AO) (IForster Schreiber et al 



20061 . 12OO9I : Ihehnert et ani2009h . now provide direct access at z ~ 2 to Ha and other rest- 



optical emission lines not strongly affected by dust, in contrast to Lya. However, optical 
spectroscopy such as that presented here has the advantage of lower sky and telescope back- 
ground emission and higher sensitivity instruments and detectors, as well as broader redshift 
coverage. 

The observations are described in § |2l evidence for kinematic information is presented in 



13. 3[ and we discuss the results in the light of current models, theory, and other observations 
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in § m Throughout the paper we adopt the currently favored cosmological parameters 
= 0.3, Vl/y = 0.7, and h = Ho/72 km s'^Mpc-^ 



Observations and Data Reduction 



2.1. Lyman Break Galaxy Sample 



Our sample was drawn from 46 candidate and confirm ed LBGs in the Hubble Deep 
Field North (HDF-N). We used the same sample selection as iLowenthal et al I (119971 ): blue 
in "5 - /" (^450 - JgM < 1-22)1 but extremely red in "f/ - B" {Usoo - ^450 > 1-41), the 
colors expected for a blue, star-forming galaxy at 2 < z < 3.5 with a Lym an continuum 
break. This color criterion is similar to that adopted by ISteidel et al. I (|2003f). b ut is slightly 
less restrictive in the blue U — B, red B — I region, which ILowenthal et al I (119971 ) found 
allowed detection of "Lya-break galaxies" at higher redshift in addition to pure Lyman 
break galaxies. Seven of the 46 sources were included via that relaxed criterion. The galaxy 
magnitudes fell in the range 24.01 < /814 < 26.27, wit h a median va l ue of 25.34. The sample 
is thus more than 1 magnitude fainter tha n that of iPettini et al I (120011 ) and more than 2 
magnitudes fainter than that of lErb et al I (120061 ). We chose both targets with confirmed 
spectroscopic redshifts and those without, since the success rate for LBGs selected in this 
way is typically so high (> 90% for bright sources, R < 25.5). The m edian reported redshift 



was z = 2.7. Many of the galaxies discussed here were studied by ISteidel et al. I (Il996bl ) 



and ILowenthal et al I (119971 ) but with lower spectral resolution and/or less attention to the 
position angle of the spectrograph slit. 

We took full advantage of the fine pixel scale and high resolution of the HDF WFPC2 
images, 0.04"pixel~^ with FWHM ~ 0'.'14. We preferentially chose LBGs that showed clear 
signs of one or more of the following: (1) significant spatial extent > 1" in the HDF-N 
image; (2) multiple knots and/or close pairings of LBGs in the HDF-N; and (3) previous 
spectroscopic evidence for some kinematic features. We also prioritized bright targets over 
faint ones. In all we observed 32 kinematic LBG targets in the HDF-N. Of the 32 total LBG 
targets observed, 21 showed some promise in the data of allowing dynamical mass estimates; 
the remaining discussion will focus exclusively on those 21 objects plus one serendipitously 
discovered LBG, which together we call the "high-priority targets". A complete list of all 22 
targets is given in Table [21 and the 11 unused targets are listed for completeness in Table [31 



^Note that we are using AB magnitudes throughout, transformed from ST magnitudes as detailed in the 
HDF information posted on the World Wide Web: to transform from ST magnitudes to AB magnitudes, 
add 1.31, 0.399, -0.199, and -0.819 to U3Q0, B45Q, Veoe, and Isi4, respectively. 
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Most of those 11 were not useful for kinematic measurements because they had insufficient 
signal-to-noise ratio, no strong emission line, or no detectable spatial extent in our spectra. 

We compared our target list to the list of LBGs h osting AGN in the HDF according to 
Chandra X-ray flux as reported by iLaird et al I (120061 ): no sources match within 10" down 
to 0.5-2 keV flux levels as low as 3 x 10~^''erg s~^cm~^ , so we conclude that none of our 
targets hosts a luminous AGN. 



2.2. Observations 

Our observations were made during six observing runs from May 1997 to April 1999 
at the 10-meter Keck-II telescope using the Low Resolution Imaging Spectrograph (LRIS; 
Oke 1995). We designed six slit masks around the LBG targets; high-priority targets were 
observed with multiple masks to increase integration time. Each target was examined closely 
in the HDF image, and then its corresponding slit was tilted to match the long axis of any 
extended emission from the target galaxy. Slits were also tilted to cover multiple targets 
in close pairs or groupings. The masks were designed so that the position angle (PA) of 
most slits was within 20° of the PA of the mask, meaning the slits were aligned close to 
perpendicular to the dispersion direction. But some slits were tilted as much as 40° from 
the mask PA and the normal to the dispersion direction. Each mask had between 23 and 
25 slits total (not including holes for alignment stars). About half the slits on each mask 
were designed around targets in the HDF and Flanking Fields for other programs. Slit 
lengths for our LBG targets ranged from 8-35", with a typical length of 15". All slits were 
I'.'l wide, although the slit tilts cause that width to be projected to a dimension smaller on 
the sky by a factor of cos(i), where i is the relative angle between the slit and the mask. 
The mask PAs were optimized for the anticipated mean hour angle of the HDF at the time 
of observation, in the sense that we tried to keep the majority of slits aligned with the 
parallactic angle to minimize slit losses due to differential refraction. We relied on our own 
astrometric measurements of the HDF and Flanking Fields to derive target positions. 

For deriving wavelength solutions, HgNeArKr arc lamps were observed through each 
mask. Strong night sky lines were used to fine-tune the wavelength solution of some spectra. 

We used a 600 1/mm grating blazed at 5000 A, resulting in a dispersion of 1.28 A pixel"^, 
a typical resolution of 4 pixels ~ 5 A full-width-at-half-maximum (FWHM) or 300 km s~^. 
This represented a compromise between a lower-resolution grating that might have delivered 
higher signal-to-noise per pixel and a higher-resolution grating that would have provided 
more detailed kinematic information. The grating angle was adjusted so that central wave- 
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lengths for slits at the mask centers were between 5000 and 6OOOA, giving a spectral coverage 
of 2620A and spectra ranging from a minimum of 3456A to a maximum of 6978A. The spa- 
tial pixel scale for slits perpendicular to the dispersion direction was 0.215 arcsec pixel"^, 
although this scale is compressed (more arcsec pixel"^) by l/(cos i) for tilted slits, where i 
is the relative PA between the slit and the mask. 

Individual exposure times ranged from 1800 sec to 2700 sec, with at least two exposures 
through each mask to help reject cosmic rays. Spatial dithering was not possible due to the 
tilted slits. Targets received between two and 17 exposures each, depending on priority and 
repeat placement on multiple masks. This resulted in total exposure times ranging from 
9,900 to 50,000 sec per target. 

The weather was mostly clear and photometric but some exposures were affected by 
cirrus. The seeing was typically 1". 

Images of all the targets are shown in Fig. [21 with the location and PA of each slitlet 
superposed. Targets observed at multiple PAs are shown once for each PA. 

We also i ncluded in our analysis some data from an earlier Keck/LRIS observing run 



in April 1996 (iLowenthal et al 1119971 ). The slits for that run's slitmasks were tilted only to 
accommodate multiple objects, not to align to the major axis of extended objects. Never- 
theless, some of those slits did fall at the same PA as the slits in the observing runs aimed 
at kinematic study, allowing us to add the datasets together. 

The new observations are summarized in Tabled! and properties of the targets are listed 
in Table H 



2.3. Data reduction 

We reduced the Keck/LRIS spectra using a combination of IRAF tasks and custom 
IRAF scripts and C routines. Raw images were bias-subtracted, cosmic-ray cleaned, flat- 
field corrected, and geometrically corrected to account for optical distortion. 

The next stage was carried out in one of two ways: (1) The wavelength solutions for 
all the slits in a single exposure were measured simultaneously, the image was rebinned so 
that all slits had the same pixel-wavelength mapping and orientation perpendicular to the 
dispersion, and then background sky emission was subtracted and individual slits extracted 
as single two-dimensional spectra; or (2) individual slits were extracted from each exposure 
and then treated as long-slit spectra: a wavelength solution was derived from a corresponding 
arc lamp image, the data were rectified, the wavelength solution was fine-tuned using night- 
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sky emission lines, and background sky emission was subtracted. 

Finally, for each slit all individual exposures were combined using exposure-time weight- 
ing and spatial registration to make a single two-dimensional, sky-subtracted spectrum. 
One-dimensional (ID) spectra of most targets were extracted from the two-dimensional (2D) 
image. Representative 2D and all the ID spectra are shown in Fig. [3] and Fig. HJ respectively. 



3. Redshifts and LBG Kinematics 

In this section we present several new redshifts in the HDF, and then describe our 
measurements of kinematic features in the reduced spectra and how we use them to derive 
estimates of dynamical masses, along with some significant caveats on those estimates. 



3.1. Redshifts 



We examined each one-dimensional spectrum visually to measure or confirm the target's 
redshift. For difficult cases we tried cross-correlating the one-dimensional sp ectrum with var- 



i ous te mplate spectra, including an average of confirmed LBG spectra from iLowenthal et al 



(119971 ) and a hi gh-signal-to- noise ratio (S/N) spect rum of the gravitationally-lensed LBG 



MS 1512-CB-58 flYee et allll996l : IPettini et al 1 120001 ) kindly provided by M. Pettini. In no 



case did the cross-correlation produce a reliable new redshift. In all, we found three tentative 
and one robust new redshifts; confirmed, refined, or constrained four more; and were unable 
to obtain useful redshift constraints for two. We were also able to measure the redshift for a 
new LBG the image of which fell on one of the slits but that was not in our original target 
list, for a total of five new redshifts. 

Table m gives the redshift obtained for every object in the sample. The redshift quality 
Qz is given on a scale from 1 to 4 such that = 1 means there is little hope of assigning 
a redshift given the signal-to-noise ratio in our data, 2 means real features are evident but 
the redshift is not secure, 3 means the redshift is probable, and = 4 means the redshift 
is definitely secure, with multiple spectral features identified. 



3.2. Kinematic Measurements 



We next searched for evidence in our LBG spectra that would allow us to constrain 
their dynamical state or even estimate their masses or their host halo masses. Measuring 
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robust dynamical masses, however, is not straightforward, and the mass estimates we derive 
are subject to several major caveats. We discuss some of the most significant of these below. 



3.2.1. Non-virialized systems: implications from simulations 

The velocities and velocity limits we measure may simply not reflect orbital motion under 
gravity. Peculiar motion of galaxies or subgalactic clumps in unvirialized systems such as 
loose groups or very young clusters could play a significant role. Chance near superposition 
of physically unrelated sources could also lead to erroneous mass estimates. Since many of 
our targets are indeed clumps of knots or close pairs, these factors may strongly affect our 
results and interpretation. 

To try to assess the significance of the non-virialization scenario, we t urn to the "phe- 



nomenological" n-body and semi-analytic ACDM galaxy formation models of lvan Kampen fc Crawford 



( 120071 ). Their "Model 3" represents a combination of two scenarios: star formation domi- 
nated by quiescent star formation in disks, usually just the central galaxy, and star formation 
dominated by merger-induced starbursts (both galaxy-galaxy and halo-halo mergers). One 
run of that model produced 1 3,257 simulated LBG s at z ~ 3 "detected" in 1 deg^ with 



criteria similar to those used by lSteidel et al. I (120031 ). We searched the catalog for projected 



close pairs of LBGs and examined their characteristics. We find 1206 pairs with separation 
r < 7". In Fig. [5l we plot the relative velocity Af in each pair vs. the projected separation 
Ar on the sky, where the velocities are calculated using the two galaxies' redshifts and their 
peculiar velocities. We focus on the simulated pairs with redshift difference A2; < 0.001 
(where "redshift" means cosmological redshift only, without considering peculiar velocity), 
which we take to be truly physically related. There are 103 such pairs in the simulated 
catalog, or 8.5% of all pairs with r < 7"; these are represented as open triangles in Fig. [51 
We also plot the values for the four close pairs in our sample (see § H]) . 

We calculate for each simulated close pair of LBGs a dynamical mass 



Mdyn = rdynt^rot/G (1) 

where frot is the maximum observed circular velocity, which in this case we take as Aw, 
and Tdyn is the radius at which the velocity is measured, in this case the separation between 
the pair. We also calculate the total (dark matter) halo mass of each pair as the sum of the 
two individual halo masses listed for the galaxies in the simulation catalog. 

Fig. [HI shows calculated dynamical mass vs. total halo mass for the close pairs in the 
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simulation of Ivan Kampen &: Crawford I (l2007l ). Surprisingly, there is no strong correlation 
seen. This implies that the velocities in fact do not (or not always) reflect true virial velocities. 
They may instead reflect infall velocities of subclumps or galaxies streaming into a potential 
well; such infall velocities are likel y larger than virial velocities by at least a factor of two. We 
conclude that, if the simulation of Ivan Kampen fc Crawford I (120071 ) is correct, the pair-wise 
velocities in our LBG sample may provide at best weak constraints on the true total masses 
of our LBG systems. 



3.2.2. Incomplete spatial sampling 



A second caveat is that we may be sampling the kinematics of only a small region, either 
central to a larger potential well or embedded but not centered in a large massive proto- 
galactic cloud or even a large but UV-faint galaxy. As discussed above in § [H several recent 
groups have used optical and NIR photometry to study the rest-optical S EDs of LBGs in 



the HDF-N, and found that their stellar masses were typically M < 10 Mrr^ (iPapovich et al 



2001 



iDickinson et a. 



Shapley et al 



2005 



2004 



Huang et al 



Sawicki fc Yee Ill998l : iRigopoulou et al 112009 : iLabbe et al II2005I : 



20051 ) . These stellar masses provide a firm lower limit to 



the total dynamical mass. 

In the local universe, rotation curves of mass ive disk galaxies rapi dly reach maximum 
rotational velocity fmax, typically within a few kpc (ISofue fc Rubin 1120011 ) . similar to the half- 
light radii of the LBGs studied here. Therefore, if LBG kinematics behave at all like those of 
local disks then we might not in fact be missing higher rotation velocities from faint extended 
regions below our detection limit (although we could be missing faint emission at the same 
velocity from larger radii, which would raise the mass estimates). Similarly, the central 
velocity dispersions a of massive local elliptical galaxies are generally excellent indicators 



of th e total galactic virial mass as measured by independent methods (jPeZeeuw fc Franx 



199ll ). Again this implies that even if LBGs reside in larger systems, the total mass may be 
well-sampled by the LBG kinematics. 

Finally, we can compare LBGs at 2; ~ 3 to luminous compact blue galaxies (LCBGs) 
at redshifts z < 1, which have many properties similar to those of LBGs such as small half- 
light sizes ri/2 < 3kpc, high luminosities L ~ L* (where L* is the characteristic luminosity 
in a Schechter luminosity function of galaxies today), blue optical colors, diverse, irregular 
morphologies, asymmetries A ~ 0.3, star formation rates 10 < SFR < 100 Mrr, yr~^, and 



narrow optical emission lines, a < 100 km s ^ ( 


Guzman et al 


Pisano et al 


2001; 


Guzman et al 


2003; 


Lowenthal et al 


2005). 
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Barton et al I (120061 ) found from deep optical imaging that more than half of the 27 
compact narrow emission line galaxies (CNELGs, close cousins of LCBGs) at 0.1 < z < 0.7 in 
their sample had sizes consistent with small local dwarfs, arguing against massive underlying 
host galaxies. Similar results derive from studies of nearby compact UV-lumi nous galaxies 
using GALEX jBasu-Zvch et al ll2007l : lHeckman et allboOsl i lHoopes et all2007h . The overall 
picture is that LCBGs at intermediate redshift are galaxies with luminosity L ~ L* but 
masses only 1/1 0th that of a typical L* galaxy, and dynamical/stellar mass ratios < 3. 

Constraints on stellar mass do not address the question of extended massive dark matter 
or gas halos surrounding the LBGs, of course. However, like meas urements of local galaxy 
halo masses using satellite galaxy kinematics (jZaritsky et al Ill993l ). our constraints on the 
dynamical masses of close pair LBG systems are insensitive to mass-to-light-ratio or stellar 
population uncertainties. 



3.2.3. Lya linewidth complications 



The third significant caveat regarding our mass estimates is that the Lya emission 
linewidths we use to estimate dynamical masses of individual LBGs and sub-clumps are 
subject to the influence of bulk inflows and outflows of gas, strong absorption by dust, and 
multiple scattering by neutral hydrogen, all of which ca n change the l ine profile. Such effects 



have been studied theoretically by severa l groups (e.g.. IWolfe 



1986; Verhamme et al 2006 



Neufeld I Il990l : iTenorio-Tagle et al 1 12004 ) . Observational evidence includes the systematic 



asymmetry and redshifting of the Lya emission line - presumably due to a c ombination of 
outflows and absorption by foreground dust - in LBGs and local starb ursts (IPettini et al 



20011 : lErb et al II2004J : iForster Schreiber et al II2006I : lErb et al Il2006f ) , and lOstlin et al I (1200 
find that the Lya emission of local starburst galaxies is spatially more extended than the 
continuum emission, presumably due to resonant scattering. We searched for but could not 
find any published or unpublished direct comparison between observed linewidths of Lya and 
those of Ha or other non-resonance nebular emission lines. Therefore our mass estimates 
based on Lya emission hne widths must be regarded as tentative. 



3.2.4. Unknown inclinations 

Finally, we generally lack constraints on the inclination angle i of our targets, so each 
mass has associated with it an unknown sin i correction factor. Apart from estimating or 
measuring inclinations of individual LBGs (as done by, e.g., Bouche et al 2007), this problem 
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can be addressed only by assuming an average inclination or else assembling a large enough 
sample of kinematic measurements that all inclinations are statistically well represented. 



3.3. Mass Estimates 

Given the caveats discussed above, we will proceed with caution in analyzing LBG 
dynamics from our sample. Here we make an empirical distinction between individual Lyman 
break galaxies, with any multiple knots of emission connected by continuous flux visible in 
the HDF images and usually separated by ~1" or less, and LBG systems, consisting of two 
or more LBGs without obvious extended emission bridging the gap(s) between them. Both 
categories may correspond to star-forming galaxies or sub-galactic clumps embedded in low- 
mass or massive dark matter halos, and both may yield useful constraints on dynamical 
mass. The LBG systems offer an additional opportunity to measure the total host halo mass 
as opposed to masses of subclu mps or individual LBGs, al t houg h given the results of our 



examination of the simulation of Ivan Kampen fc Crawford I (120071 ). we will not assume that 
such results are robust. 

We examined visually the two-dimensional spectra of all 22 of the high-priority targets 
for evidence of any of the following: (1) spatially extended emission lines, especially Lya 
but also including C IV and He II; (2) spectrally unresolved, resolved, or multiple-peaked 
emission lines, especially from LBGs with multiple knot morphologies; (3) extended contin- 
uum against which extended absorption lines might be detected; and (4) emission and/or 
absorption lines from close pairs of galaxies. The features detected are listed in Table HI 

For each LBG or LBG system with evidence of kinematic information, we then estimated 
the dynamical mass. For simple gravitational rotation, the dynamical mass would be M^yn 
as in Eq. [1] above. 

We divide our mass estimates into two categories: (1) masses estimated from extended 
emission or absorption or multiple sources; and (2) masses estimated solely from the Lya 
emission line velocity width and observed half-light radius. The measured velocity and size 
parameters and derived mass estimates for our sample are listed in Table [51 

We compare in Fig. O the distribution of velocities, which range from < 60 to nearly 
5000 km s^^, and projected separations o f the four close pair systems with the prediction 



of Model 3 of Ivan Kampen &: Crawford I (120071 ). The pair with the largest Av appears to 



lie outside the locus of simulated pairs with small Az (triangle symbols; assumed to be 
physically associated, rather than just chance projection), but there are some simulated 
pairs with small Az and even larger Av than that target - i.e., two LBGs within a few 
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xlOO kpc of each other but with very high relative velocity, Af ~ 10^ km s~^. To assess the 
probability Ptruc that each pair in our sample is a true physical association, we can compare 
the density of small-A^; pairs (triangles) to the density of large-Az pairs (dots) in different 
sections of the simulated Ar — Ai; plane. We divide Fig. into eight equal segments bounded 
by 1, 4, and 7" and 0, 2000, 4000, 6000, and 8000 km/s. Pairs with separation Ar < fare 
excluded because we would be unlikely to resolve such pairs spatially given the typical seeing 
in our data. We calculate the probability Ptrue = -^true/^^totai, where A^tme and iVtotai are the 
number of small-A^; and total pairs, respectively, in each section. We find that Ptrue ranges 
from 57% for the two real LBG pairs (x symbol) in the lower left corner of Fig. [5] to 6% for 
the pair in the upper middle section. In other words, according to the simulation, two of the 
four pairs have a better than 50% chance of being true physical associations, while one has 
a 94% probability of being a chance superposition. The value of Ptrue for each pair is shown 
in Fig. El 

For each of the eight LBGs with Lya emission (six of which are also included in the 
spatially extended emission category), we fit a gaussian function to the emission line and 
measure the half-light radius ri/2 in the HDF J814 image. We do not deconvolve the half-light 
radii with the WFPC2 PSF FWHM =0'.'14. This will have only a small effect on the larger 
sources, but for the smaller sources our measurement will overestimate the true half-light 
radius and therefore the enclosed dynamical mass. We did smooth the HDF images with 
a gaussian filter in an attempt to constrain (and correct for) the true seeing by matching 
the seeing-blurred spatial FWHM in the 2D spectra. However, significant uncertainties in 
the measured FWHMs (especially for low S/N sources) and even variation within single 2D 
spectra made this approach impractical. 

We deconvolve each observed linewidth with the measured 5A spectral resolution to 
derive intrinsic linewidths. Based on the discussion above of dynamical studies of local and 
high-redshift galaxies, we adopt the FWHM of the line profile as a reasonable but untested 
proxy (modulo any outflow and extinction effects) for circular velocity frot, we adopt ri/2 
for Tdyn, and we calculate an estimate, again based on Eq. [H of the dynamical mass within 
the half-light radius. These measurements and the linewidth-based masses are also listed in 
Table 0, and all the derived masses are plotted vs. redshift in Fig. [71 Six of the eight Lya 
emission lines were at best marginally resolved, so we plot the derived dynamical masses for 
those sources as upper limits. 



3.4. Notes on Individual Objects 

Here we discuss the results for each target, listed in order of increasing RA, as in Table [21 
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3.4- 1- hd4-025 9^1947 
This relatively bright, elong ated source has a di sky isophote and captivated our interest 



as a potential high-redshift disk. I Cohen et al I (120001 ) report for their object No. 2 a redshift 



z = 0.904, with low confidence {Quality = 9, where 11 is the lowest confidence). Our 
spectrum shows no sign of Mg II absorption at that redshift in either 2D or ID, and no 
C III (although at 3634A, our spectrum has low sensitivity). We conclude that the redshift 
assignment z = 0.904 seems unlikely. 

However, the final coadded spectrum shows weak emission lines at 4499, 5869, and 
6045 A, plausibly matching [O II]3727, H/3, and [O III]5007, respectively, at redshift z = 
0.207. We therefore conclude (though with low confidence, given the lines' weakness) that 
the source is a low-redshift interloper that slipped through our color selection filter. 



3.4.2. hd4_l 076^1 847 



This compact, hi gh-surface brightness source is the brightest target in our survey. 
Lowenthal et al I (119971 ) were unable to m easure a redshift, bu t found weak evidence for Mg II 
absorption at z = 1.0155 and z = 0.879. I Cohen et al I (I2OOOI ) rep ort z = 0.882 with lo w con- 
fidence {Quality = 9, where 11 is the lowest confidence), and cite lCohen et al I (119961) as the 



redshi ft source, but we are unable to find the target in that reference. iFernandez-Soto et al 



(I2OOII ) argued that the published spectrum does not support that re dshift, and added that 
their photometric redshift technique favors 2;phot = 0.00. Meanwhile, iBudavari et al I (120001 ) 
report a photometric redshift using HST/WFPC2+NICM0S ^phot = 2.67. 

We accumulated 24.5 ks of integration time observing the source with LRIS. Never- 
theless, we find no strong emission or absorption lines, and we are unable to obtain an 
unambiguous redshift for the source, nor can we measure any kinematic signature. 

No strong continuum break is visible in our spectrum, which extends down to 3320A 
with continuum detectable by eye in the 2D image down to 3920A. Since no strong Lya 
emission or absorption is seen redward of that wavelength, we conclude that the emission 
redshift must be z < 2.17. 

We again find weak evidence for Mg II absorption at z = 1.0155 and z = 0.879, 
as reported in iLowenthal et al I (119971 ). Unfortunately the Mg II lines for the lower red- 
shift coincide with night sky lines, complicating their detecti on and measureme nt. The 
lower redshift, which is consistent with the z = 0.882 value of I Cohen et al I (l2000l ). is sup- 
ported by possible Fe II absorption. If real, either or both sets of redshifted absorption 
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lines could be due either to intrinsic absorption from the emitting galaxy itself, or to inter- 
vening absorption if the emitting galaxy is at higher redshift. The redshift of the galaxy 
G OODS J12364 85+6 21203.4, with a position centroid 2'.'2 from hd4_1076_1847, is reported 



bv ICohen et al I mm to be z = 1.010. GOODS J123640.85+621203.4 has a highly linear 
morphology extending at least I'.'S in the HDF, aligned to within 10° of the separation vector 
between the two sources. The close agreement between its emission redshift and the tenta- 
tive Mg II absorption redshift we find in the spectrum of hd4_1076_1847 suggests that it 
is either responsible for the absorption, or else associated with another as yet unidentified 
absorbing galaxy - perhaps hd4_1076_1847 itself. 

We are therefore unable to confirm any of the previously published redshifts or derive 
a new one for this source, apart from the lower and upper limits mentioned above. No 
kinematic information is available from our spectrum. 



3.4.3. C4-09 



Object C4-09 in the catalog of ISteidel et al. I (Il996bl ) is a remarkable source consisting 
of four bright k nots of emis s ion w ith nearly identical colors, all within an area barely more 
than 1" across. IZepf et al. I (119971 ) investi gated and fina l ly reie cted the source as a possible 
gravitational lens system. As reported by lSteidel et al. I (jl996bl ). the spectrum shows strong 
emission with two peaks at 5118.7 and 5132. 6A. Continuum is clearly detected redward of 
the emission but not blueward, supporting the interpretation of the emission lines as Lya 
at redshifts 3 . 211 and 3.222, respectively, slightly lower than the z = 3.226 reported by 



Steidel et al. I (1l996bl ). Weak absorption lines are seen at 5314, 5489, 5624, and 5633A, 
corresponding to Si II, O I, C II, and Si IV, respectively at an average absorption redshift 
z = 3.216, between the redshifts of the emission line peaks. 

We observed C4-09 with two PAs, one (98°) aligned with the longest axis of the paral- 
lelogram of knots, and the other (137°) closer to the short axis. The four knots visible in the 
HDF image are so close together, with spacings between 0'.'57 and 1'.'17, that the continuum 
and the red side of the emission line appear unresolved in both of our two-dimensional spec- 
tra. The blue side of the Lya emission line, however, appears to be slightly extended spatially 
in the long-axis spectrum, with FWHM=1'.'7, compared with the seeing FWHM=l'.'l, and 
offset towards the east with respect to the redder line and the continuum. No such spatial 
extent or offset is visible in the short-axis spectrum. No sign of ordered rotation is visible. 



The spatial extent and offset of the blue peak with respect to the red peak and the 
continuum are consistent with a scenario in which some subset - one, two, or three knots - 
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of the four emission knots is responsible for the blue peak, and another subset at a different 
velocity gives rise to the red peak. The blue and red Lja peaks are well fit in the spectral 
direction by gaussian profiles with FWHM 5.5 and 7.0 A, respectively, i.e., unresolved or 
marginally resolved. After deconvolution with the 5 A resolution, those widths correspond 
to < 135km s~^ and 290km s~^, respectively. The 14 A separation between the two peaks, 
however, corresponds to a velocity difference of 820 km s~^. If the wavelength separation is 
indeed due to radial velocity differences among two or more of the four knots (rather than, 
e.g., an absorption line superposed on a broad emission line), then we can use it to estimate a 
system dynamical mass. Given the maximum spatial separation (which will yield maximum 
mass) of 1'.'17, we derive a dynamical mass from Eq. [DMdyn < 1-3 x lO^^M©. 

Alternately, the double-peaked emission profile is also well-fit by a single narrow absorp- 
tion line with A = 5125. 7A and FWHM = SA superposed on a broad emission line with 
A = 5126. 5A and FWHM = lAA. Such a broad emission line could be caused by strong 
gas outflows, and the absorption line could be caused by a foreground clump or screen of 
gas, even dust-free, that resonantly scatters Lya photons out of the line center. In this 
scenario, the emission profile FHWM of 800 km s~^ would reflect outflow velocity, perhaps 
coupled with resonant scattering, rather than circular velocity of a virialized system. The 
redshifts from the Lya emission and absorption lines would then he z = 3.2164 and 3.2170, 
in excellent agreement with the redshift z = 3.216 from the inters tellar absorption line s but 



contrary to the commonly red-shifted Lya emission reported by IShapley et al I (120031 ) and 



Erb et al I (120041 ). which they cite as evidence of outflows in LBGs. Thus an outflow model 



for C4-09 appears somewhat problematic. 

We are not able with the current data to distinguish definitively between the two sce- 
narios. If outflow is indeed responsible for the line profile, then the true dynamical mass is 
most likely significantly lower than the value derived above. 



3.4.4. hd4_0818A037 



This source is compact and relatively bright, but continues to elude efforts to obtain a 
definitive redshift. 



Cohen et al I ( 1l996l ) list a redshift z = 2.268 and I Cohen et al I (l2000l ) report a cor- 
rected redshift z = 2.500 with high confidence {Quality = 2, where l=highest confi- 
dence). However, we see no Al II1671 A absorption matching either redshift, only possible 
O 11302 A/Si II1304A a nd possible N V1243 A in P-Cygni profile matching z = 2.500; no 
C IV1550 A is detected. iBudavari et al I (120001 ) give a photometric redshift 2;phot = 1-44. 
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Continuum is clearly detected in our spectrum with a spatial profile matching an unre- 
solved PSF. No strong emission or absorption lines are detected. Weak emission lines may 
be present at 5925 and 5976 A. If these were [O III] at 2; = 0.194, we would expect HP at 
5804 A and/or [O II] at 4450 A, but neither is seen. Possible weak absorption lines appear 
at 4566, 4678, 5069, and 5447 A, but none is well-detected, nor were we able to discern 
a pattern among them indic ating a source redshift. We find no support for the tentative 



redshift z = 2.04 reported by iLowenthal et al I (119971 ) based on cross-correlation, apart from 



the possible line at 5069 A, which would correpond to Al II1671 A. 

Without solid absorption or emission line features available, we also find no kinematic 
indicators in our spectrum. 



3.4.5. hd4_0298_0744 



hd4_0298_0744 is another elusive target . ICohen et al I (120001 ) report z = 2.801 with 
high confidence, while iBudavari et al I (I2OOOI ) find Zphot = 1.76. Continuum is clearly visible 
in the two-dimensional spectrum even before sky subtraction down to the blue limit of our 
spectrum, with no strong break visible. No strong emission or absorption features are seen in 
the final two- or one-dim ensional extr a cted s pectrum. We were therefore unable to confirm 
the redshift reported by I Cohen et al I (120001 ). which we would have expected to produce a 
strong Lya emission or absorption line or continuum break visible at 462lA. Because no 
strong continuum break is visible down to 3900A in our spectrum, we conclude that the 
redshift is most likely z < 2.2. 



We cannot extract any dynamical information on the source. 



3.4.6. C4-O6 

C4-06 from the catalog of Steidel et al. (1996) is one of the brightest, largest sources 
i n our sample. The HDF image shows a linear structure - dubbed the "Hot Dog Galaxy" 



(IBunker II2OOII ) - roughly 1 x 3" in extent, with two major knots and several sub-knots. The 



LRIS mask slit was aligned along that structure. 

Strong, spatially resolved continuum is clearly detected in the two-dimensional spec- 
trum, although it appears as a single spatially continuous source, rather than the separate 
knots seen with HST resolution. The spatial extent of the continuum emission, measured 
by compressing the two-dimensional spectrum into a one-dimensional spatial profile, is 3'.'2, 
consistent with the HDF image. 
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No emission lines are visible, but several strong absorption lines are detected matching 
redshifted Lya, Si II, O I, C II, and Si IV. The Lya line is broad (78A FWHM), while the 
other lines are unresolved. Additional marginal detections include Si II, N I, and Si III. The 
redshift derived from averaging the strong narrow absorption line redshifts is z = 2.794, while 
cross-correlating the spectrum with that of a template aver a ge of 1 2 LBGs yields z = 2.802, 



similar to the redshift z = 2.803 reported by ISteidel et al. I (Il996bl ). 



The narrow absorption lines show no sign of any velocity shift across the extended 
continuum (see Figure [3]). To quantify any subtle velocity gradient, we extracted three one- 
dimensional spectra at the middle and the two ends of the linear source and intercompared 
them. The central extraction aperture was 12 pixels (2'.'6) wide and the apertures at the two 
ends were 6 pixels (1'.'3) wide, centered 8 pixels apart, so that the two end apertures had no 
overlap. 

To constrain the velocity shift between the two end apertures, we selected regions of 
the spectra containing strong absorption lines and free from residual noise from strong sky 
emission lines and cross-correlated them. The cross-correlation yielded a strong peak at SX = 
O.SA. Assuming a minimum relative wavelength accuracy of I.OA, the velocity shift between 
the two apertures is then 6v < 60 km s~^, consistent with zero within our measurement 
errors. The dynamical mass limit implied by the observed constraint on line-of-sight velocity 
gradient 6v < 60 km s~^ over 3'.'22 is Mdyn < 2.6 x lO^M©, the lowest dynamical mass 
estimate in our sample (modulo, of course, the unknown inclination and three-dimensional 
morphology of the source). 

We interpret the lack of observed velocity gradient over more than 20 kpc as evidence 
that we are observing not an edge-on disk with ordered rotation, but rather a truly linear, 
perhaps filamentary source. The object may of course be collapsing or even expanding per- 
pendicular to the line of sight; collapse is more natural under most current galaxy formation 
scenarios. 



3.4.7. hd4_l 006^0680 



This galaxy appears as a very close pair in the HDF image, with a separation of only 
0'.'3. I Cohen et al I (I2OOOI ) report a redshift for the source z = 2.969. Our mask slitlet was 
aligned along the same PA, 73°, as the two emission knots; however, 1" ground-based seeing 
certainly blurs the two into a single unresolved source in our LRIS observations. 

The spectrum shows very strong emission at 4823.2 and weaker emission at 4812. 3A. 
The stronger line appears extended both spatially and spectrally: the spatial FWHM of the 
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line is 1'.'2, comparable to the seeing, but the emission profile's faint wings extend to cover 
3'.'3. The spectral FWHM is 7.0A, compared to the resolution of SA, and the total detected 
emission spans 26. 4A, including the two peaks and their wings. 

The weaker line's spatial position centroid is offset 0'.'25 from that of the brighter line, 
closely matching the separation of the pair in the HDF image. 

Weak continuum is detected both redward and blueward of the emission lines, but 
appears to be stronger on the red side. No other emission or absorption features are seen. 

It is highly unlikely that this double emission line corresponds to an optical line emitted 
by a low-redshift {z < 1) system. The [O II] doublet would, at z = 0.293, be separated by 
3.5A, not the llA we observe. Any [O III] or Balmer line should be accompanied by other 
optical lines as well. We conclude that the most likely interp retation is Lyo: at z = 2.967, 



confirming with only a slight revision the redshift reported by I Cohen et al I (120001 ). 



The double emission line profile may be due to a single emission line with strong ab- 
sorption superposed by an intervening cloud or galaxy. Alternatively, as for C4-09, it may be 
caused by radial velocity differences between the two emitting knots. Given the close match 
between the spatial separations in the image and in the spectrum, we adopt the latter sce- 
nario. The velocity difference between the emission peaks then corresponds to 677 km s~^, 
while the entire 26.4 A span of the emission line corresponds to ~ 1650km s~^, and the spec- 
tral FWHM of 7.0 A corresponds, after deconvolving with the SA instrumental resolution, 
to 4.9 A, or ~ 300 km s'^ 

To derive a mass estimate for the system, we assume each knot in the pair emits one 
emission line, and we therefore adopt r = 0'.'3 (2.3 kpc) and a = 677 km s~^ to obtain 
M = 2.4 X 10^^ M0. The large total range in observed emission velocity, 1650 km s~^, may 
be due to bulk gas outflow induced by merging of two or more sub-units and/or by the 
subsequent starburst. 



3.4.8. hd4_l 486^0880 



Despite the relative brightness of this target, its redshift remains elusive. iLowenthal et al 

(119971) repo r ted a tentative redshift z = 2.47 based on a few possible a bsorption lines, 
honipson I (|2003| ) estimated a photometric redshift Zphot = 2.80, while iBudavari et al 



( 120001 ) give none. The source is smooth, slightly extended, and elongated in the HDF image. 



Continuum is easily detected along the entire length of our spectrum. However, no emission 
lines are visible, nor are any strong absorption features. There are possible absorption lines 
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at 4386, 4528, 4794, 4816, 5168, and 5573A, but we were unable to discern any convincing 
pattern matching our redshifted template spectra or our previous tentative redshift z = 2.47. 
Cross-correlation with the template spectra likewise revealed no robust redshift. The lack of 
a strong continuum break redwar d of 4000A i r aplies z < 2.29, contrary to the photometric 
redshift 2;phot = 2.80 calculated by [Thompson I (120031 ). 

With no strong absorption or emission features, we were also unable to derive any 
kinematic information from the spectrum. 



3.4.9. hd4_0367_0266 

The HDF image shows the source hd4_0367_0266 to be an elongated object with two 
bright knots separated by 0'.'6 and an extended tail terminating in a fainter knot, all nearly 
aligned. The redshift of ^ = 2.931 reported by Lowenthal et al. (1997) was based on several 
absorption lines. The total extent of the source shown in the HDF is 2'.'3; the LRIS mask 
slitlet was aligned to cover all three knots as well as the tail. 

Our spectrum shows weak extended continuum emission with FWHM ~ 3'.'4. No emis- 
sion lines are detected. Several possible weak absorption lines appear in the extracted one- 
dimensional spectrum at 4956, 5117, and 5478A. These lines would correspond to Si II, O I, 
and Si IV, respectively, at the previously reported redshift z = 2.931. A possible broad 
absorption trough at 4757A could be a damped Lya absorber at z = 2.913. Overall our 
confidence in the redshift is somewhat lower than reported in Lowenthal et al. (1997), and 
we downgrade the redshift quality to = 3. 

To search for velocity gradients or shifts across the galaxy, we extracted two spatially- 
independent one-dimensional spectra centered 2'.'4 apart and cross-correlated them, using 
only the clean regions of the spectra that were free from sky noise residuals and showed evi- 
dence of absorption features. The cross-correlation function displays a weak peak consistent 
with 6v = 0, but unfortunately the signal-to-noise ratio was too low to allow any robust 
measurement of kinematics. 



3.4.10. hd2_1928A041 (C2-06) and C2-05 



The source hd2_1928_1041 (C2-06 in the catalog of Steidel et al 1996) is separated by 
only 2'.'0 spatially from C2-05, and we consider them here as a close pair. Each of the galaxies 
is a comp act source with one br ight knot and a small cloud of extended emission as seen in 
the HDF. Isteidel et al. I Jl996bh reports z = 2.845 for C2-05, but this is revised to z = 2.005 
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bvlCohen et all fl2000h . For hd2_1928_1041, IPapovich et al I (1200 If ) report z = 2.009, citing 



Lowenthal et al I (119971 ): we here slightly revise that redshift. 



Both targets are covered in a single LRIS slitlet on our slit masks. Strong, unresolved 
continuum is detected from both objects. No emission lines are visible, but numerous absorp- 
tion lines corresponding to interstellar species such as Si II, O I, C II, Si IV, C IV, and Al II 
are well detected in the two spectra, corresponding to redshifts z = 2.005 (hd2_1928_1041) 
and z = 2.008 (C2-05). The velocity difference implied by the redshifts is 6v = 390 km s^^. 
Cross-correlation of absorption line regions in the two one-dimensional extracted spectra pro- 
duces a strong peak at Sv = 364 ± 87 km s~^ (where the error is estimated by dividing the 
FWHM of the cross-correlation peak by 10), consistent with the simple redshift difference. 

The observed velocity difference and spatial separation provide a dynamical mass esti- 
mate Mdyn ~ 5.8 X IO^^Mq for the two-source system. 



3.4.11. hd2_1739A258 



A highly elongated, multiple-knot source embedded in extended emission, extending 
over at least 2'.'6 along the long dimension, to which our LRIS slitlet was aligned. 



Spectroscopic redshifts reported in the litera ture include z, 



1 -pec—. V 

et al I I2OOOI ) with high confidence 



2.72 ( Lowenthal et al 



19971 ) with low confidence a nd Zspen = 1.980 (| Cohen 
{Qual ity = 3, where l=best). iFernandez-Soto et al I ( 



mm cite z. 



spec 



2.002, and Cowie et al 



( 12004) cites z = 1-98, but we could not find the source of either measurement in the litera- 



ture. 



the source. 



Wirth et al. I (120041 ) reported that a redshift was not measurable in their spectrum of 



^Photometric red shifts calculated fo r hd2_1739_1258 include ^^phni^ = 1.640 (|Fernandez-Soto et al 

19991 ). 2phot = 2.854 (iThompson II2003I ). and ^phot = 1.34,2.07 (iBudavari et allboool ). 



Strong, extended continuum emission is detected in our 39 ks integration with a spatial 
FWHM ~ 9 pixels, or 1'.'9 and a total extent of 3'.'9. The continuum dies below about 
4075A, implying z ~ 2.35 if the dip is due to Lyo; blanketing. However, we also detect two 
possible weak absorption lines at 5305, 5318A, which match Mg II at z = 0.897. These 
absorption hues could be due to intervening gas, but they could also be intrinsic to the 
target. We dete ct no features that sup port either the redshift z = 1.980 (ICohen et al II2OOOI ) 
or Zspcc = 2.72 (ILowenthal et al 1119971 ). Given the lines we do detect weakly, we tentatively 



conclude that hd2_1739_1258 in fact lies at 2; = 0.897 rather than at z > 2. 



We extracted two spatially independent 1-dimensional spectra from the extended two- 
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dimensional spectrum and cross-correlated them in an attempt to measure velocity gradients 
across the source, but no statistically significant cross-correlation signal was detected. 



3.4.12. hd2A398A164 
This gal a xy co nsists of a single compact core with two short wings of extended emission. 



Cohen et al I (1200(1 ) report z = 2.237 



Our LRIS slitlet was aligned with the long axis of the system. The spectrum shows a 
single weak emission line at 3934. 2A with FWHM ~ 6 A, i.e. barely resolved, with clear 
continuum redward but none blueward. Interpreting the line as Lya, which is supported 
by the continuum break, we derive a redshift z = 2.236. No other strong emission lines are 
detected, although possible weak C IV and Al II absorption and He II emission are visible. 
The width of the line corresponds to a deconvolved velocity width of 240 km s~^, which we 
adopt as an upper limit. 

Neither the continuum nor the strong emission line shows any spatial extent in our 
two-dimensional spectrum, so we are unable to constrain kinematics for this system. 



3.4.13. C3-02 



phot 



1.720, and Budavari et al 



J2000h 



C3-02 is a small source with a single compact bright knot of emi ssion and a dif f use tai l 
extending barely 1", along which our LR IS slit was roughly ori e nted. ISteidel et al. I (Il996bl ) 
reported a tentative redsh ift z = 2.775. IPernandez-Soto et al ( 1999 ) report a photometric 
redshift z- 



report Zphot = 2.218. 



Spatially unresolved continuum is easily detected with average S/N~ 11 per resolution 
element over the entire length of our 25.5 ks spectrum, which extends from 3650-6285A. 
No emission lines are visible. Several possible weak absorption lines appear at 4317, 4663, 
and 4608A, the first of those being the strongest (i^VFobs ~ SA). None of those wavelengths 
matches features expected from galaxies aX z = 2.775. Both the photometric and spectro- 
scopic redshifts previously reported therefore seem unlikely to be correct. There is a possible 
match to O 11302, Si II1304, C II1335/1336, and Si IV1402/1403 at z = 2.316. No Lya 
emission is visible at the expected wavelength for that redshift, but the continuum shape is 
consistent with a moderate break there. Given the weakness of the lines, we adopt this red- 
shift only tentatively [zq = 2). Kinematic measurements of this source are also impossible 
from our spectrum. 
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3.4- U- hd2_0698A297 

hd2_0698_1297 consists of two distinct knots embedded in an asymmetrical blob of diffuse 
emission. The knots are O'.'S apart, and the diffuse emission visible in the HDF F814W image 
extends for 1'.'2. The system is 2'.'7 from hd2_0705_1366, an unresolved point source (measured 
FWHM =0'/14) in the HDF image. 

Because of slit mask orientation constraints, we observed this pair at two different 
position angles: for 25.5 ks at 42°, and for 10.2 ks at 58°. Due to the close proximity of the 
two sources, the I'.'l-wide slitlets at both PAs include emission from both targets. 

Both spectra show weak continuum with S/N ^ 1 per resolution element and a single 
emission line from each source. The emission line from hd2_0698_1297 is stronger in the 
spectrum with PA=42°, while that from hd2_0705_1366 is stronger at PA=58°, presumably 
because of slit placement with respect to each source. We measure the wavelengths of the 
emission lines from the stronger spectrum in each case: 5396A (hd2_0698_1297) and 5312 
(hd2_0705_1366), each with FWHM = 8A, corresponding for red shifted Lya to z = 3.43 9 



and 3.370, respectively, within 6z < 0.01 of the values reported in iLowenthal et al I (119971 ). 
There is evidence for an absorption trough bluewards of the emission line in hd2_0698_1297, 
supporting our identification of the line as Lya. We cannot confidently rule out the possibility 
that each emission line is in fact spectrally unresolved, so we adopt the measured FWHM 
as an upper limit in each case. 

The continuum and the emission lines are all spatially unresolved, except perhaps con- 
tinuum from hd2_0698_1297, which shows a slight spatial extent FWHM ~ 1'.'5, although 
this may be due at least partly to small spatial registration errors in co-adding the 18 indi- 
vidual exposures. 

Given the wavelength difference 84 A and projected separation 2'.'7 between the two 
sources, and assuming the wavelength difference to reflect gravitationally-induced dynamics, 
we derive a system dynamical mass Mdyn = 1 x IO^^Mq, the largest estimated dynamical 
mass in our sample. 



3.4.15. hd2_0705A366 



See hd2_0698_1297 (above). 
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3.4-16. hd2_0529A567 

See hd2_0624_1688 below. 

3.4.17. hd2_0624J688 

This triple-knot compact source is part of a close pair with hd2_0529_1567, 6'.'25 away. 
Both were covered by the shtlet at PA = 28.9° that target e d hd2 _0725_1818 and hd2_0743_1844 



see below). hd2_0624_1688, for which iLowenthal et al I (119971 ) reported z = 2.419, was in 



our original kinematic sample target list, while hd2_0529_1567 was not: at Isu^ab = 27.233, 
it offered little hope of providing a useful spectrum. 

However, our 2D spectrum shows clear continuum and a single strong emission line from 
each source. The continuum emission from hd2_0624_1688 shows a significant drop blueward 
of the emission line. Continuum from the fainter hd2_0529_1567 is of course much weaker, 
and the proximity of the emission line to the blue end of the spectrum prevents us from 
assessing any continuum break across the line. 

The single emission line from hd2_0624_1688 lies at 4162. 7A, yielding r edshift z = 2.424 



assuming that the emission line is Lya, compared with z = 2.419 reported by ILowenthal et al 



( 119971 ). No other absorption or emission lines are visible, and Lya is the only plausible in- 
terpretation. The emission linewidth is FWHM = 6.3 A, corresponding to 3.8 A or 
~ 275 km s^^ after deconvolution with the 5 A resolution, and we adopt that value as an 
upper limit. 

The emission line from hd2_0529_1567 lies at 3973. 3A, yielding z = 2.268, again under 
the assumption that the line is redshifted Lja. The deconvolved linewidth is FWHM = 
3.5 Aor ~ 270 km s~^, which we adopt as an upper limit. No other absorption or emission 
features are seen. 

Together with the projected separation of 6'.'25 (50 kpc), the observed velocity difference 
of 1394 km s~^ yields a dynamical mass estimate for the pair system Mdyn = 2.2 x lO^'^M©. 



3.4. 18. hd2_0725J818 



This small (ri/2 = 0'.'26) source lies only 1'.'3 from hd2_0743_1844. ILowenthal et al 



(I1997I ) reported a spectroscopic redshift z = 2.233 for hd2_0725_1818 based on a Lya emission 



line, but only a tentative redshift z = 2.39 for hd2_0743_1844. The two sources have similar 
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colors, with hd2_0743_1844 bluer by only 0.27 mag in B^^o — 1^14 and redd er by only 0.05 



mag i n [/300 — -B450, less than the spread in the colors of LBGs reported by IPapovich et al 



(120051 ): hd2_0743_1844 is fainter by 0.8 mag in J814 (all measurements from the catalog of 



WiUiams et al 1996). 

Our LRIS spectrum covers not only both of those sources, b ut als o a spiral galaxy to 



the north with a redshift z = 1.148 reported by iPhillips et al I (119971 ) and, to the south, 
hd2_0624_1688 and hd2_0529_1567 (see above). 

Continuum from both galaxies is easily visible in the two-dimensional spectrum, but 
with some overlap of their spatial profiles, given their mere 1'.'3 separation. The flux from 
hd2_0725_1818 is 2-3 times stronger than that hd2_0743_1844, consistent with their pho- 
tometric measurements. Several absorption features are apparent in the two-dimensional 
image; while they seem to span spatially the continuum from both sources, the weakness of 
the fainter source makes this appearance difficult to confirm visually. 

We extracted a one-dimensional combined spectrum of the pair. To assess their redshifts 
and velocities independently, we also extracted two separate one-dimensional spectra in 
the following way: First we extracted a one-dimensional spectrum of the brighter source, 
hd2_0725_1818, using only those pixels in the two-dimensional spectrum extending spatially 
from the peak of the source away from hd2_0743_1844, i.e. the side of the spatial profile that 
is less contaminated by fiux from the fainter source. We then used this ID extraction to 
model the entire 2D spectrum and subtracted the scaled model from the 2D image. This left 
a clean, isolated spectrum of hd2_0743_1844, which we then extracted to ID. We repeated the 
process in reverse, modelling and subtracting the 2D spectrum of hd2_0743_1844 to produce 
a clean, isolated and uncontaminated spectrum of hd2_0725_1818. 

The absorption lines in the ID spectrum of hd2_0725_1818 are well matched by inter- 
stellar Sill, CII, and Alll, and s tellar SilV, CIV a t z = 2.232, thus confirming and slightly 



revising the redshift reported by iLowenthal et al I (119971 ). We also find a tentative detection 



of C III in emission, similar to their earlier result. Lya lies at the blue end of the spectrum, 
which is clean but low S/N; no emission or break is visible in either 2D or ID. 

The spectrum of the fainter source, hd2_0743_1844, also shows absorption corresponding 
to Si II, C IV, and O I. No Al II is detected, though the S/N is low in that part of the 
spectrum. The implied redshift is z = 2.230 ± 0.004, consistent with that of hd2_0725_1818, 
and we adopt that redshift with quality = 3 rather than 4 given the source's extreme 
faintness. 



To measure any small velocity difference between the two sources, we cross-correlated 
the independent ID spectra against each other, using only those parts of the spectra free 
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of bright sky emission line residuals. The result was a strong correlation peak (correlation 
strength = 0.2) at 6v = 9 km s~^, consistent with 6v = km s~^ within our nominal velocity 
resolution of 60 km s~^. For comparison, we also cross-correlated each ID spectrum with a 
ID spectrum of blank sky extracted from the same 2D image. No strong correlation peak is 
seen anywhere near 6v = 0, with a maximum of only 0.017 within 3000 km s~^ of Sv = 0, 
and a maximum correlation strength 0.1 for the whole spectrum. We thus conclude that 
the correlation between hd2_0725_1818 and hd2_0743_1844 is significant, and that they have 
identical redshifts within our uncertainties. 



Adopting an upper limit 6v < 60 km s ^ and the projected separation I'.'S or 10.9 h ^ 

f 

'■djn 



kpc, we then derive a dynamical mass upper limit Mdyn < 8.6 x IO^Mq, the lowest derived 



mass in our sample. 

3.4.19. hd2_0743J8U 

See hd2_0725_1818 above. 

3.4.20. hd3A633A909 

An elongated, wispy source consisting of two clumpy wings extending over 2'.'1, with 
a O'.'S gap separating them, and aligned roughly east- west. The source's total is o phota l 



colors in the HDF Vesion 1 catalog satisfy the 5-dropout criteria of iLowenthal et al I (119971 ). 
However, in the HDF Version 2 catalog, the two sides' colors are significantly different, with 
the western side bluer by 1.30 mag in -8450,^5 — Veoe.AB and 0.52 mag in Veoe.AB — -^814,ab- 
Furthermore, in the Version 2 catalog photometry, neither component individually satisfies 
either the [/-dropout or the 5-dropout criterion. 

Our LRIS slit covered both wings, as well as the source hd3 1824_1945 7'.'4 away and 



two other sources, one on either end of the slit. I Cohen et al I (120001 ) report z = 4.050 for the 
source. 

Only extremely very faint continuum is detected at the position of hd3_1633_1909 in the 
2D spectrum, with no absorption lines discernible. A single emission line appears, however, 
at 6I45A. The emission line, though apparently spectrally resolved with FWHM = 8.4A, 
is not well fit by the [O II] doublet at z = 0.649; the peaks of the two components would be 
separated by 4.1A, while the line in the spectrum is well fit by a single gaussian. If the line is 
Lya, the source's redshif t is z = 4.056 , the h ighest redshift in our sample and consistent with 



the redshift reported by I Cohen et al I (120001 ). The resolution-deconvolved velocity width of 
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the emission line is FWHM = 330km s ^. 

No spatially extended kinematic information is available from our spectrum. 



3.4.21. hd3A 824A945 

This compact dou ble- knot sour c e was covered by the same slit that included hd3_1633_1909 
[see above), 7'.'4 away. I Cohen et al I (120001 ) report Zsp = 2.050 (object 145 in their online list). 



Our spectrum shows strong continuum with spatial FWHM =V!6 extending without 
any obvious break virtually to the blue limit of the spectrum at 3750A. No strong emission 
lines are visible; only a few weak possible absorption lines are present. We see no compelling 
evidence to support the redshift z = 2.050 reported by I Cohen et al I (I2OOOI ): for example, 
Al II1671A, C IV, and Si IV, which often appear in strong absorption in LBG spectra, are 
not seen. 

The lack of a strong Lya continuum break down to 3 75 A implies z < 2.2, ir i con- 
fiict with the photom etric redshifts Zphot = 2.30 calculated by iFernandez-Soto et al I (I2OOII ) 
and 2;phot = 2.395 by iBudavari et al I (120001 ). Meanwhile, the lack of strong emission lines 



including [O II] 3727A implies z > 0.75. 

No kinematic information is available from our spectrum. 



3.5. Summary of results 

Of the 22 targets in our high-priority sample, we confirm that 14 are LBGs with z > 
2 and = 3 or 4, of which 13 provided some evidence of kinematics that we use to 
estimate dynamical mass or upper limits thereon. Two candidates are probably interlopers 
at ^ < 1, and the remaining six show insufficient features for confident redshift (or kinematic) 
measurement. The failure rate, 6/22 = 27% (or 8/22=36% if we include the two interlopers), 
is higher than the typical failure rate for LBGs with R < 25.5, < 10%. This must be due to 



our color selection function, which is somewhat different from that of ISteidel et al. I (120031 ) 
and may be less restrictive, as mentioned above; our magnitude limit, which is 1-2 mags 
fainter than most previous spectroscopic surveys of LBGs; our use of the 600 1/mm grating, 
which provides less spectral coverage than the lower-resolution gratings commonly used for 
LBG redshift surveys; or a combination of those effects. 



Eleven of the 13 LBGs with kinematic signatures show spatially extended emission or 
absorption: eight are in close pairs, two contain three or four knots of emission each, and one 
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is a clumpy elongated source extending over 3"or ~ 25 kpc. A different subset comprising 
eight of the 13 are sources with Lya emission lines whose linewidths may provide some 
constraint on dynamical mass. The redshifts of the combined subsets range from z = 2.005 
to z = 4.056 with a median of z = 2.424. The angular extent of emission ranges from 
0.12-6"(~ 1 — 50h~^ kpc), and the adopted circular velocities range from < 60 km s~^ to 
4700 km s^^, with a median of 330 km s~^. 

We divide the mass estimates into two overlapping categories based on the source 
morphologies and angular extent: individual LBGs, comprising 11 sources, and LBG sys- 
tems, comprising the four close pairs. We derive dynamical masses for the LBGs (using 
Lya linewidths and extended and multiple-knot sources) ranging from < 4.3 x lO^M© to 
1.3 X IO^^Mq, with a median M^yn = 4.2 x l O^'^MfD. We see no evidence in any of our tar- 



gets fo r ordered rotation such as measured by iForster Schreiber et al I (120061 ) and lErb et al 
(120061 ) for some of their more luminous star- forming galaxies at lower redshift {z ^ 2). Dy- 
namical mass estimates for the LBG systems range from < lO^^M© to ~ lO^^Mo, with a 
median mass M = 1O^^M0. We emphasize that all these values are subject to the caveats 
discussed in § 13.21 and that two of the four close pairs have probability > 70% of being 
chance pr ojection rather than physical ass ociations, according to our comparison with the 
models of Ivan Kampen fc Crawford I (120071 ). 



The tentative dynamical masses we derive for individual LBGs, < AxlO^h ^MfD to 1.1 x 
lO^^/i~^M0, are very similar to those measured for brighter targets by lForster Schreiber et al 



((20061) based on Ha kinematics measured with VLT/SINFONI, mdyn ~ 0.5 - 25 x 10^°Mq, 



implying that either Lya and Ha emission linewidth measu rements are both robust or else 
they both suffer from complications of comparable scale. iForster Schreiber et al I (120061 ) 
also extrapolate the velocities beyond the half-light radii to estimate the viri al mass of a 
typical halo, finding mhaio ~ IO^^'^^'^'^'^Mq, similar to the halo masses inferred by lAdelberger 



(j2005al ) based on clustering. In contrast, the median dynamical mass 1.1 x lO^^M© that we 
measure for our four close-pair systems is an order of magnitude larger, and the maximum 
sy stem dynamical mass estimate is > lO^'^/i-^M©, more than 5a above the mean halo mass 
of iForster Schreiber et al I (120061 ). We note, however, that in no case in our sample is the 
derived dynamical mass less than the reported stellar mass for the same galaxy (§ 14.21) . which 
is reassuring. 



4. Discussion 



Having measured redshifts and searched our spectra for dynamical signatures in 13 LBGs 
and LBG systems, we now compare those results to several other observational results and 
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to theoretical expectations in order to explore implications for galaxy formation scenarios. 



4.1. Mergers, close pairs, and implications for galaxy formation 



Many lines of evidence point to the importance of merging in galaxy evolution over the 
entire cosmic history of galaxies. In this section we explore the possibility that our sample 
of LBGs includes a significant merger fraction, and we compare the characteristics of our 
sample directly to the predictions of galaxy formation theories, including masses and close 
pair fractions. 

Cold accretion of gas and dark matter not associated with major majors is also in- 
creasingly appreciated as a significant source o f fresh material for star formation and mass 
buildup (T lCook et al II2009I : lAgertz et al II2009I ). But major mergers, where the components 
have roughly equal mass, can erase disk morphologies and produce starbursts, AGNs, and 
elliptical gala xies, and even minor mergers can induce bursts of star formation or d isturbed 
morphologies (IBarnes fc Hernquist lll992l : lMihos &: Hernquist lll996l : lLin et al 1120071 ). Recent 
works have measured the evolution of the rate of galaxy mergers, the fract ion of close pairs 



of ga l axies, and the dist r ibution of galaxy morphologies over cosmic time ( Kampczvk et al 



2007 



2006a 



Bridge et al 



2007 



Lin et al II2004I : iBell et al II2006I : iKartaltepe et al 



2007 



Lotz et al 



Conselice et alll2003l ). Essentially all of those results are consistent with hierarchical 



galaxy formation models, in which galaxy mass is assembled largely through merging. 

A s noted above , early studies of LBGs recognized th e high rate of close pairs and multiple 



knots (ICollev et allll996l . 119971 : iLowenthal et allll997h . Of the 14 confirmed LBGs in our 



sample, 12 are in close pairs or have multiple-knot or highly extended, clumpy morphologies. 
All 12 are therefore excellent candidates for merging systems. Indeed, every close pair 
in our sample yielded potential kinematic information; conversely, the targets that failed 
to provide such information were predominantly low-surface brightness, elongated, diffuse 
sources, which may themselves be interacting or merging systems. Our entire sample, of 
course, is strongly biased away from isolated, symmetrical morphologies by virtue of our 
visual selection process, so it is to be expected that the sample would show a higher fraction 
of merging systems than do magnitude-limited LBG samples. 

Our sample contains only four close pair systems, so it is subject both to small number 
statistics and especially to the assumption that the pairwise velocities indicate dynamical 
mass (see § 13.21) . The components of the four pairs may be subclumps in various stages of 
infall and virialization, in which case their relative velocities could be dominated by random 
rather than circular motions. They may, however, also reflect a true, large range in the 
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masses of the LBG host halos. 

The theoretical imphcations of observed LBG close pair incidence and velocity dis- 



(Somerville et al 


2001 


Wechsler et al 


J2OOI 


) 



Wechsler et al 2001; Zhao et al 2002, 2003 



et alll200lh . 



Wechsler et al rj2ooih found that different models predicted very different halo occupation 



distributions, i.e. the number of LBGs per dark matter halo. In most ACDM models, the 
number of LBGs predicted per halo increases with mass, with typically ~ 1 LBG halo~^ for 
M ~ IO^^Mq and R < 25.5. The number of close pairs of galaxies - defined as galaxies 
within a given angular separation within a redshift interval Az < 0.04 - is an especially 
sensitive discriminator between models. Specifically, compared to the observed number of 
close LBG pairs, models postulating a one-to-one correspondence between LBGs and mas- 
sive dark matter halos ("massive halo" models) underpredict close pairs by 1.5cr, while the 
"colliding halo" model overpredicts close pairs by 4a within 15". The best agreement was 
provided by their "collisional starburst models", in which LBGs host both ongoing quies- 
cent star formation and rapid bursts of star formation triggered by minor and major galaxy 
mergers. 



Wechsler et al 


( 


2001 


), we 


satisfy the 


Wechsler et al 



criterion Az < 0.04. Fig. 6 in that paper shows the 
fraction iVpairs/iVgaiaxies as a function of pair separation for different models. All four of the 
pairs in our sample have separation < 15", so they all fall in the first bin in the figure. 
If we assume that we have observed all the close pairs in the HDF survey volume in that 
separation bin, i.e. the total number of close pairs with Az < 0.04 is iVpairs = 2, and further 



that the total number of galaxies in the survey volume with comparable selection criteria 
is roughly the number of LBG candidates in the HDF times the success rate we observe 
here, or Ai'gaiaxies ~ 46 x 64% = 29, then we find A^pairs /^galaxies ~ 0.07. This value is 
already marginally inconsistent with (higher than) the values, all ~ 0.04, predicted by the 
co llisional starburst, con stant efficiency quiescent, and accelerated quiescent models in Fig. 6 
of lWechsler et al I (1200 ll ): the massive halo model underpredicts close pairs by a factor of more 
than lOcT (including only the error bars of Wechsler et al 2001). The colliding halo model, 
on the other hand, overpredicts it by a factor of less than la. Obviously, with such small 
numbers of close pairs and large uncertainties on the selection function and completeness of 
our sample, the observed incidence ratios should be viewed as only tentative. 

The close pairs and multiple knot systems we study here all have maximum angular 
separation r < 7", and there f ore a ll fall in the first bin in the separation-pair fraction 
plot (Fig. 6) of I Wechsler et al I (I2OOII ). where indeed the discrimination among models is the 
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strongest. Because of the biased nature of our target selection, we do not attempt to compare 
quantitatively the observed pair fraction with the predicted one. But we do note that close 
pairs are represented over the entire range of dynamical masses we derive, 10^° — IO^^Mq. This 
distribution is likely due either to (1) misinterpretation of the kinematic signatures (i.e. the 
velocities are not due primarily to orbital motion under gravity), resulting in erroneous mass 
estimates or (2) misidentification of individual knots of emission as separate galaxies, thus 
complicating the comparison with the models; or else (3) it is inconsistent with the general 
prediction of I Wechsler et al I (120011 ) that the number of LBGs per halo should increase with 
halo mass. 



Shu et al I (1200 ll ) use observed sizes and SFRs of LBGs to conclude from their models 
that LBGs must be hosted by halos with circular velocities Vc = 100 — 300 km s~^, although 
for some combinations of gas-to-total and dust-to-total mass ratios the distribution can reach 
400 km s~^. If LBG close-pair dynamics are indeed indicative of host halo mass, then the 
range of halo masses in our sample is sig ni ficantly larger and the typical mass significantly 
higher than predicted by IShu et al I (120011 ) . IZhao et al I (120021 ) further explore the possibility 
of using the observed pairwise velocity distribution (PVD) to discriminate between different 
galaxy formation models, but their investigation focuses on scales between 100 kpc and 
30 Mpc, much larger than the 10-kpc scale of multiple-knot and close pair targets studied 
here. 



Zhao et al I (l2003al ) analyze dark matter halo A^-body simulations and report that mass 



accretion occurs in a rapid phase, during which the halo's circular velocity is also predicted to 
rise rapidly, followed by a slow phase, during which the circular velocity remains almost con- 
stant. However, they examine only isolated halos with total mass fixed at M = 10^^h~^MQ, 
for which the maximum circular velocity predicted is ~ 300 km s~^ and the maximum halo 
radius is ~ 30h^^ kpc. These spatial scales are roughly comparable to those of three systems 
in our sample: two close pairs and the highly extended source C4-06. The velocity scale is 
comparable in only one close pair system. Most of the LBGs and LBG systems in our sample 
are therefore not apparently well described by the models of IZhao et al I (l2003a( ). However, if 
we interpret our results within their framework, the large range in system dynamical masses 
we derive could imply that we are seeing the systems at different evolutionary stages, some 
before the rapid mass accretion phase is completed and the circular velocity plateaus, and 
some after. Typical group-sized halos today have circular velocities ~ 500 km s~^, and clus- 
ter halos have circular velocities ~ 1000 km s~^. Only one of our LBG system velocities lies 
well above that range, with Av = 4700 km s~^; the rest have velocities consistent with those 
of normal galaxies to normal clusters of galaxies today. As noted above in § 13.31 (Fig. 
the distribution of pairwise velocities ar id projected separations of our fo ur pairs is consis- 
tent with the predictions of Model 3 of Ivan Kampen &: Crawford I (120071 ). supporting their 
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scenario of quiescent disk star formation combined with merger-induced starbursts. 



The dynamical masses we derive span the range of masses of dwarf galaxies to giant 
elliptical galaxies today. Our modest sample is certainly biased and is not representative of 
all LBGs at 2 ~ 3. The range of estimated dynamical ra a sses i s similar to that derived for 
simulated LBGs in Model 3 of Ivan Kampen fc Crawford I (120071 ) (§ 13.31) . but the kinematics 
of both simulated and real datasets may in fact not reflect true virial masses, especially 
given the lack of correlation with combined halo mass in the simulation (Fig. E]). However, 
taken together, the sizes, SFRs, morphologies, velocities, and derived individual and system 
dynamical masses of the LBGs in our sample imply a large range of physical scenarios and 
dark matter halo masses. They also further suggest that, rather than being identified one- 
to-one with the centers of future massive spheroids, bulges, or elliptical galaxies, LBGs are 
merging sub-clumps within halos that are destined to become galaxies with a large range of 
masses and potentially of morphological type as well. This scenario is consi stent with hybrid 



galaxy forr nation models such as the favored collisional starburst model of ISomerville et al 
( I2OOII ) and IWechsler et al I (120011 ). but not with their massive halo model. 



4.2. Kinematic vs. Stellar Masses 

Here we compare our kinematic mass estimates to the stellar masses of the same galaxies 
calculated by SED fitting. 



Papovich et al I (120051 ) expanded on the sample of iPapovich et al I ( l200ll ). fitting stellar 
population models to SEDs of UV-bright galaxies (i.e. LBGs) at 1.9 < 2 < 3 in the HDF 
based on fluxes from HST/WFPC2 and NICMOS. The stellar masses for the 53 galaxies with 
S/N > 20 in the NICMOS F160W band range from 1 x 10^ - 3 x W^Mq, with a median 
around 3 x lO^M©. The stellar mass estimates for 765 of 1682 NICMOS-detected sources with 
spectroscopic or photometric redshifts have since been supplemented by photometry from 
Spitzer/IRAC and were kindly provided to us by C. Papovich (2007, private communication). 

We searched the 5j>zfeer-supplemented NICM0S-I-WFP2 catalog for positional matches 
to our high-priority LBG sample. All 22 targets have at least one match well-detected by 
NICMOS within 0'.'6. We visually examined and confirmed the position of each match in the 
HDF F814W images. 

Of the 22 LBG targets, 20 have stellar mass estimates in the Papovich catalog, the other 
two lacking reliable redshifts. Neither of those two, however, is among the 13 LBGs in our 
sample showing kinematic signatures. Therefore we can directly compare all the individual 
and system (i.e., halo) dynamical masses we estimated to the stellar masses estimated by 
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Papovich et al I (120051 ) . The HDF Nicmos Mosaic (HNM) catalog numbers, redshifts used, 



and stellar masses of all matched objects are listed in Table [61 

Six targets have two matches in the Papovich catalog within 1'.'3. Of those, three 
are in our list of 13 with dynamical mass estimates: C4-09, C4-06, and hd3_1633_1909, each 
producing two matches. Both matches in each case correspond to emission k nots we consider 



to be part of the same LBG. For C4-09 and C4-06, the redshifts used by iPapovich et al 



(120051 ) to calculate SED fits and stellar masses were f or both matche s essen tially the same 



redshifts we report here. For hd3_1633_1909, however, iPapovich et al I (l2005l ) used the same 
redshift we find {z = 4.056) for only one of the two matches, while for the other they used 
z = 2.175, i.e. what we consider to be a single LBG they consider to be two sources at widely 
disparate redshifts. The resulting stellar mass estimate for the low-redshift source (HNM 
1499) is ~ 1/3 that of the high- redshift source (HNM 1506). We consider below the stellar 
mass estimate of only the high-redshift source. The low-redshift source certainly contributes 
flux to our spectrum. However, as our dynamical mass estimate of hd3_1633_1909 is based 
only on the Lya emission line width, the effect on our mass estimate should be negligible. 
On the other hand, the redshift z = 2.175 of the low-redshift source may be in error; if the 
correct redshift is in fact close to 2; = 4.05, then the stellar mass estimate for hd3_1633_1909 
would increase by roughly a factor of two. 

We make the following three comparisons of the stellar mass estimates to our dynamical 
mass estimates: (1) for the eight Lya linewidth-based dynamical mass estimates or upper 
limits, we compare those values directly with the single matching stellar mass estimates 
(seven sources) or the two components of a double match (one source, C4-09) for a total of 
nine comparisons; (2) for the three individual but extended or clumpy sources with single 
or multiple matches in the NICMOS catalog, we combine the stellar mass estimates and 
compare the total to our individual LBG dynamical mass estimated from extended absorption 
or emission velocity or, if neither of those is available, from Lya linewidth; and (3) for the 
close pairs, we combine the two stellar mass estimates and compare the sum to the system 
dynamical (halo) mass from the pair's relative velocity. The four-knot source C4-09, which 
yields both Lya linewidths and a double emission line peak, presumably from different 
emission knots, has two NICMOS matches, so we apply both comparison types 1 and 2 (i.e. 
we separately compare the smaller dynamical and stellar masses - perhaps due to individual 
knots - and the combined stellar and system dynamical masses). These comparisons are 
shown in Fig. [HI 

The individual LBG stellar masses (comparison types 1 and 2 above) range from 2.50 x 
10^ to 1.6 X lO^°/i-^M0 with a median of 1.6 x IO^Mq. Two sources show stellar masses 
roughly equal to their dynamical masses. The stellar mass range and median, however, are 
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both less than the range and median of dyn amical masses; the ty pical mass ratio M^yn/M^ ~ 
5, comparable to those of LCBGs at z < 1 (IGuzman et al II2003I ). The stellar masses a re als o 
small compared to those measured for 181 LBGs in the EGS by iRigopoulou et al I (120061 ). 



10^ ^^Mq with median 2.95 x 10^°Mq. No obvious corr elation is seen iii Fig. E ] between our 



dynamical mass estimates and the stellar masses from Papovich et al ( 20051 ). although 
trend could be masked by the six upper limits on dynamical mass from Lya linewidths. 



Erb et al I (120061 ) likewise compared stellar and dynamical masses for 68 galaxies at 
lower redshift z ^ 2, finding a weak correlation between the two, with a mean stellar mass 

= 3.7 X IO^^Mq - more than 20 times the median stellar mass for our fainter, higher- 
redshift sample from iPapovich et al I (120051 ) - and a typical mass ratio Mdyn/M* ~ 2 vs. 5 
for our sample. They also found that ~ 15% of their sample had large dynamical-to-stellar 
mass ratios M^yn/M^ > 10 (vs. 10/12 = 83% for our sample), suggesting that this was hkely 
due to young galaxy ages that had not yet built up significant stellar mass, or else faint older 
stellar populations that escaped detection. 

The same scenarios may explain the large dynamical-to-stellar mass ratios and lack 
of correlation w e observe. Our sample is more than one magnitude fainter than that of 
Erb et al I (120061 ). so we are probing farther down the LBG luminosity function. Local galax- 
ies and even starburst galaxies show a wide range in mass-to-light ratio, from M/L ~ 0.01 
in solar units for extreme starbursts s uch as H II galaxies t o M/L ~ 10 for quiescently star- 
forming galaxies like the Milky Way (iGuzman et allll996l ). The large range is presumably 
due to a combination of galaxy age (older systems showing stronger NIR flux from stars), 
burst age (younger systems showing stronger UV flux from stars), merging history, and dust 
obscuration. The large range we observe in dynamical-to-stellar mass ratio may in fact re- 
flect a wide range of intrinsic LBG properties including dark matter halo or sub-clump mass, 
star formation history, merging history and status, and burst age. 

We also compared the system masses derived from the four close pair relative velocities 
with the combined stellar mass estimates from Papovich et al (2005; comparison type 3 
above). The right-hand panel of Fig. [8] shows the comparison. Here the range in combined 
stellar mass is smaller than for the individual LBGs, while the range in system dynamical 
mass is larger. The average total stellar mass is 5.0 x IO^Mq, and the average system 
dynamical-to-stellar mass ratio is 1.1 x lO'', ranging from Mdyn/M* < 3 to Mdyn/M^, = 
2.6 X 10^. This large mass ratio range suggests that one or both of two scenarios are at play: 

(1) the LBGs are embedded in dark matter halos with a wide range of masses that are not 
strongly correlated to the luminosity and star formation rates of the LBGs they host; and 

(2) the relative velocities of the close pairs may reflect various stages of infall, i.e. not all 
the close pairs are necessarily in virialized systems (see § 13.21) . 
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5. Summary and Conclusions 

We have studied candidate and confirmed LBGs at 2; ~ 3 in the HDF-N with deep 
two-dimensional spectroscopy, finding 3 tentative and 2 robust new redshifts. We also find 
evidence for kinematics in 13 LBGs and four LEG close-pair systems. The observed mor- 
phologies and dynamics, including numerous cases of close pairs and multiple sources, appear 
to be most consistent with collisional starburst models of galaxy formation, rather than mas- 
sive halo models with only one LEG per halo. We derived tentative dynamical mass estimates 
from the observed velocity signatures, finding a range from 4.3 x lO^Mo to 1.3 x lO^^M© 
with a median Mdyn = 4.2 x IO^^Mq for the individual LEGs and a range from < lO^^M© 
to ~ IO^^Mq with a median mass M = IO^^Mq for the close-pair systems. All the mass 
estimates are subject to strong caveats including unknown inclination corrections, Lya ab- 
sorption, outfiows, and non-virialized systems. All our dynamical mass estimates are larger 
than the corresponding stellar mass measurements. However, we find no evidence for a cor- 
rel ation betw e en ste llar mass and dynamical mass, in contrast to the weak correlation found 



by lErb et al I (120061 ) for their sample of brighter, lower-redshift sources. 



We argue on the basis of comparison with local compact starburst galaxies, other studies 
of distant star-forming galaxies, and galaxy formation theories and simulations that the sizes, 
SFRs, morphologies, velocities, and dynamical mass estimates imply a wide range of physical 
scenarios and mass scales giving rise to LEGs at z ^ 3. 

Given the large time investment at Keck needed to observe even this small sample, 
further progress in constraining the dynamical masses of LEGs will require technological 
advances including the deep diffraction-limited imaging and spectroscopy promised by James 
Webb Space Telescope and increasingly available through adaptive optics and NIR integral 
field spectroscopy at ground-based 10-meter class telescopes. 
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Table 1. Observing Runs at Keck II Telescope 



Date (UT) 

1997 May 5-6 

1997 May 30 - June 1 

1997 June 6 

1998 April 28 - May 1 

1998 May 23-25 

1999 April 15-18 



Note. — The Low 
Resolution Imaging 
Spectrograph with a 600 
1 mm~^ grating blazed 
at 5000A and slit widths 
of I'/l were used for all 
observing runs, giving 
spectral resolution about 
300 km s"^ FWHM (130 
km s"-*^ gaussian a). 
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Table 2. Lyman Break Galaxy Kinematic Targets in the Hubble Deep Field. 



Namo=' HDFv2 ^ RA (2000)'^ DEC (2000)'= I814.ab'^ PA" ti„t ' 

mag ° ks 



hd4_0259. 


.1947 


4-916 


12 


:36; 


:38, 


,605 


62: 


12: 


;33, 


,83 


25. 


,080 


110, 26.0 


21.6 


hd4_1076. 


.1847 


4-878 


12 


:36; 


:40, 


,963 


62: 


12: 


;05. 


,30 


24. 


,009 


26.0, 137.0 


24.5 


C4-09 




4-858 


12 


:36; 


:41 


,245 


62: 


12: 


;03. 


,07 


24. 


,904 


98.0, 136.9 


21.6 


hd4_0818. 


.1037 


4-445 


12 


:36; 


:44, 


,640 


62: 


12: 


;27. 


,39 


24. 


,586 


190.0, 26.0 


24.5 


hd4_0298. 


.0744 


4-316 


12 


:36; 


:45, 


,087 


62: 


12: 


;50. 


,81 


24. 


,983 


145.0 


24.5 


C4-06 




4-555.1 


12 


:36; 


:45, 


,409 


62: 


11: 


;53. 


18 


24. 


190 


43.0 


25.5 


hd4_1006. 


.0680 


4-289 


12 


:36; 


:46, 


,947 


62: 


12: 


;26. 


,08 


26. 


,082 


73.0 


9.9 


hd4_1486. 


.0880 


4-382 


12 


:36; 


:46, 


,951 


62: 


12: 


;05, 


,34 


24. 


,957 


13.0, 26.0 


25.5 


hd4_0367. 


.0266 


4-52 


12 


:36; 


:47, 


,720 


62: 


12: 


;55, 


,79 


25. 


,579 


26.0, 108.0 


21.6 


hd2_1928. 


.1041 


2-454 


12 


:36; 


:48, 


,266 


62: 


14: 


18. 


,42 


25. 


,081 


165.0 


25.5 


C2-05 




2-449 


12 


:36; 


:48, 


,338 


62: 


14: 


16. 


,63 


24. 


249 


165.0 


25.5 


hd2_1739. 


.1258 


2-585 


12 


:36; 


:49, 


,811 


62: 


14: 


15. 


18 


24. 


,399 


6.5, 148.0 


38.9 


hd2_1398. 


.1164 


2-525 


12 


:36; 


:50, 


,123 


62: 


14: 


;01. 


,03 


25. 


,596 


83.0 


9.9 


C3-02 




3-550 


12 


:36; 


:51, 


,335 


62: 


12: 


;27. 


,51 


25. 


,706 


174.0 


25.5 


hd2_0698. 


.1297 


2-604 


12 


:36; 


:52, 


,451 


62: 


13: 


;37. 


,84 


25. 


,787 


42.0, 58.0 


41.9 


hd2_0705. 


.1366 


2-637 


12 


:36; 


:52, 


,760 


62: 


13: 


;39. 


,09 


25. 


,916 


42.0, 58.0 


41.9 


hd2_0529. 


.1567 


2-751 


12 


:36; 


:54, 


,205 


62: 


13: 


;35. 


,83 


27. 


,233 


28.9, 31.0 


38.9 


hd2_0624. 


.1688 


2-824 


12 


:36; 


:54, 


,627 


62: 


13: 


;41. 


,37 


26. 


139 


28.9, 31.0 


38.9 


hd2_0725. 


.1818 


2-903 


12 


:36; 


:55, 


,077 


62: 


13: 


;47. 


,00 


25. 


,406 


28.9, 31.0 


50.0 


hd2_0743. 


.1844 


2-916 


12 


:36; 


:55, 


,184 


62: 


13: 


;48. 


,07 


26. 


214 


28.9, 31.0 


50.0 


hd3.1633. 


.1909 


3-853 


12 


:36; 


:58, 


,967 


62: 


12: 


;23. 


,42 


25. 


,708 


77.0 


21.6 


hd3.1824. 


.1945 


3-875 


12 


:37; 


:00, 


,127 


62: 


12: 


:25. 


22 


24. 


,671 


77.0 


21.6 



^Target names are as in lLowenthal et al ] i 19971 1 and arc in one of two forms: (1) hdn_xxxx_yyyy, 
where n is 2, 3, or 4 and represents the WFPC2 CCD chip in which the source falls, and xxxx and 
yyyy the pixel coordinates of the source on that chip in the HDF Version 1 catalog, which was 
used for original target selection; or, (2) for some of the brightest targets, Cn-Om from fs'teidel et al. | 
jl996b|) . Targets are listed in order of increasing RA. 

bource name from HDF Version 2 catalog of I Will iams et al matched to our target by 

positional proximity and visual confirmation via the HDF images. 

'^J2000 equatorial coordinates of each target from the HDF Version 2 catalogs [Williams et al | 
{T99I). 

'^/gl4 AB isophotal magnitude from the HDF Version 1 catalog, used for original target selection. 
^Position angle(s) (degrees east of north) of slitlct used to observe target. 
^Total integration time on target at all PAs. 
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Table 3. Unused Kinematic Targets in the Hubble Deep Field. 



Name-'' HDFv2 ^ RA (2000)= DEC (2000)'^ I814,ab'^ 

mag 



hd4_0460. 


.1146 


4- 


■491 


12 


:36; 


:43, 


,253 


62: 


12:38.85 


25. 


,680 


hd4_1994. 


.1406 


4- 


■631 


12 


:36; 


:45, 


,310 


62: 


11:38.47 


26. 


,269 


hd2_1949. 


.0599 


2- 


■239 


12 


:36; 


:45, 


,886 


62: 


14:12.09 


25. 


342 


hd4_1341. 


.0299 


4- 


■85 


12 


:36; 


:49, 


,515 


62: 


12:20.11 


26. 


103 


hd3.0408. 


.0684 


3- 


■ 243 


12 


:36; 


:49, 


,814 


62: 


12:48.80 


25. 


,985 


hd2_0664. 


.0879 


2- 


■373 


12 


:36; 


:50, 


,302 


62: 


13:29.73 


25. 


,951 


hd2_1410. 


.1282 


2- 


■594 


12 


:36; 


:50, 


,683 


62: 


14:03.15 


26. 


,098 


hd3.0378. 


.1536 


4- 


■ 677 


12 


:36; 


:51, 


,620 


62: 


12:17.31 


25. 


,979 


hd3.0457. 


.2023 


3- 


■915 


12 


:36; 


:53, 


,107 


62: 


12:00.74 


25. 


,770 


hd2_0434. 


.1377 


2- 


■643 


12 


:36; 


:53, 


,422 


62: 


13:29.53 


25. 


,457 


hd3.1455. 


.0430 


3- 


■ 118 


12 


:36; 


:54, 


,727 


62: 


13:14.72 


25. 


141 



^Target names arc as in iLowenthal et al | j 19971 ) and arc in one of two 
forms: (1) hdn_xxxx_yyyy, where n is 2, 3, or 4 and represents the WFPC2 
CCD chip in which the source falls, and xxxx and yyyy are the pixel coor- 
dinates of the source on that chip in the HDF Version 1 catalog, which was 
used for original target selection; or, (2) for some of the brightest targets, 
Cn-Om from fs'teidel et al. | jl996b| ) . Targets are listed in order of increasing 
RA. 

^Source name from HDF Version 2 catalog of IWilliams et aTl il996|) , 
matched to our target by positional proximity and visual confirmation via 
the HDF images. 

^J 2000 equatorial coordin ates of each target from the HDF Version 2 cat- 
alogs IWilliams et al I il996D . 

'^7gl4 AB isophotal magnitude from the HDF Version 1 catalog, used for 
original target selection. 



Table 4. Redshifts and Spectral Features 



Foaturc{s)'' 



hd4.0259.1947 


D.2D7 


2 








new z 




hd4.1076.1847 


I.DID < z < 2.17 


1 






_ 


Possible Mg II abs 




C4-09 


3.226 


4 






Lya 


quadruple knot 




hd4.0818.1037 


{2.5DD) 


1 












hd4.0298.0744 


< 2.2 


1 












C4-06 


2.794 


4 






Si II, O I, C II, Si IV 


"Hot Dog" 




hd4.lD06.0680 


2.969 


3 






Lya 


double knot and em. peak 




hd4.1486.0880 


< 2.29 


1 












hd4.0367.0266 


2.931 


3 












hd2.1928.1041 


2.005 


4 


Si 


II, 


O I, C II, Si IV, C IV, AI II 


C2-06; pair with C2-05 




02-05 


2.008 


4 


Si 


II, 


O I, C II, Si IV, C IV, AI II 


pair with hd2.1928.1041 




hd2.1 739.1258 


0.897 


2 








new z 




hd2.1398.1164 


2.236 


3 






Lyci 






C3-02 


2.316 


2 








z = 2.775 unlikely; new z 




hd2.0698.1297 


3.439 


4 






Lyci 


pair with hd2.0705.1366 




hd2.07D5.1366 


3.370 


4 






Lyci 


pair with hd2.0698.1297 




hd2.0529.1567 


2.268 


4 






LyQ 


new z; pair with hd2.0624_ 


.1688 


hd2.0624.1688 


2.424 


4 






LyQ 


pair with hd2.0529_1567 




hd2.0725.1818 


2.232 


4 






O I, Si II, AI II 


pair with hd2.0743.1844 




hd2.0743.1844 


2.230 


3 






O I, Si II, AI II 


new z; pair with hd2.0725_ 


.1818 


hd3.1633.1909 


4.056 


3 






Lya 






hd3.1824.1945 


0.75 < z < 2.2 


1 













^Targets arc listed in order of increasing RA, as in Table 

'^Redshift derived from Lya emission line wavelength or stellar and/or interstellar absorption line wavelengths (for 
2 > 1) or optical nebular emission lines (for z < 1). Typical uncertainties are 0.001 in z. Values in parenthesis indicate 
redshifts published by others that wc were unable to confirm or refute. For details on redshifts listed as ranges or 
limits, see text. 

"^Redshift quality. — no spectrum obtained; 1 —highly uncertain, no definitive features identified; 2 —real features 
are evident but the redshift is not secure; 3 — redshift probable, with two or more spectral features identified; 4 — redshift 
definite, with multiple spectral features identified. 

'^Spectroscopic feature(s) used for measuring dynamical mass. Lya is in emission only; all others are in absorption 
only. 



Table 5. Dynamical Mass Estimates 







Mass from Lycv Lincwidth 






^lass from 


Extended Emission 


Namc^ 


V ( FW H M)^ 




A / J ^ 




''dyn^ 


'■dyn'' 


A^dyn' 




(km s ^ ) 


{aresec) {/i ^ kpc) 




{km s ^ 


) (arcsec) 


(h^l kpc) 




hd4_0259.1947 












- 


- 


hd4_1076.1847 


_ 




_ 


_ 




- 


- 


C4-09 


< 135 


0.14 1.0 


< 4.3 X lO" 


820 


< 1.2 


< 8.5 


< 1.3 X 10^^ 




290 


0.14 1.0 


1.9 X 10^" 






- 


- 


hd4_0818.1037 


_ 








_ 


- 


- 


hd4_0298.0744 








_ 


_ 






C4-06 


- 


- 


- 


< 60 


3.2 




< 2 X 10^^ 


hd4_1006.0680 


300 


0.22 1.7 


3.4 X 10^° 


677 


0.3 


2.3 


2.4 X 10^^ 


hd4_1486.0880 
















hd4_0367.0266 




- 


- 








- 


hd2_1928.1041 








390 


2.0 


16.5 


5.6 X lO" 


C2-05 








390 


2.0 


16.5 


5.6 X 10^^ 


hd2_l 739.1258 
















hd2_1398.1164 


< 240 


0.26 2.2 


< 2.8 X 10^" 










C3-02 
















hd2_0698.1297 


< 350 


0.45 3.3 


< 9.0 X 10^" 


4707 


2.7 


19.6 


9.8 X 10^^ 


hd2_0705.1366 


< 350 


0.14 1.0 


< 2.8 X 10^" 


4707 


2.7 


19.8 


9.8 X 10^^ 


hd2_0529.1567 


< 270 


0.12 1.0 


< 1.6 X lOl" 


1394 


6.3 


50.7 


2.2 X 10^^ 


hd2_0624.1688 


< 275 


0.31 2.5 


< 4.2 X 10^" 


1394 


6.3 


50.1 


2.2 X 10^^ 


hd2_0725.1818 








< 60 


1.3 


10.6 


< 8.6 X lO" 


hd2_0743.1844 








< 60 


1.3 


10.6 


< 8.6 X lO*" 


hd3.1633.1909 


330 


0.64 4.5 


1.1 X 10^^ 










hd3.1824.1945 
















^Targets arc listed in order of 


increasing RA, as in Table f2l 












''FWHIVI of Ly 


a emission line 


measured from the ID spectrum and deconvolved vi'ith 


the spectral resolution FW H M 


= SA. Typical 


uncertainties are 


30 km s ^ ^ . 














■^Half-light rad: 


ius as measured 


in the HDF images. Uncertainties are typically 0. 05. 









As in c, but transformed from arcscconds to kpc using our adopted fiat cosmology with Q.r,-i — 0.3, — 0.7, Hq — 

72 h km s~ Mpc~ "'" . Uncertainties are ^ 0.4 kpc. 

^Dynamical mass calculated from A/^iyn — ''"1/2''^-?' W-ff M ^ /G. Uncertainties {based solely on uncertainties in '(^FVV/fA/ and and 
neglecting uncertainties due to non-virial motions, extended stellar or dark matter components, source inclination, outflows, etc.) on 
the values in this column are 10^ Mq. 



Observed velocity spread of kinematic feature{s). For close pairs, the velocity separation is entered twice, once for each component. 
Typical uncertainties are 30 km s""'". 

^Effective spatial size, spread, or separation as measured in the HDF images. Uncertainties are typically 0. 1. 

^ As in g, but transformed from arcseconds to kpc using our adopted flat cosmology with — 0.3, — 0.7, Hq — 

72 h km s~ "'" Mpc~ . Uncertainties are ^ 1 kpc. 

^Dynamical mass calculated from -M^y^ — r^jy^Au^/G. For close pairs, the mass represents the mass of the system and is entered twice, 
once for each component. Uncertainties (based solely on uncertainties in Aij and i^dyn and neglecting uncertainties due to non-virial 
motions, extended stellar or dark matter components, source inclination, outflows, etc.) on the values in this column are ^ 2 X 10^. 
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Table 6. Stellar Masses from Papovich et al (2005) 



Namc^ 


HNM id'' 














(arcscc) 




(Mq) 


hd4.0259.1947 


748 


0, 


,01 


-1.000 


0.00 


hd4.1076.1847 


1078 


0, 


,02 


0.882 


2.20E9 


C4-09 


1114 


0, 


,57 


3.216 


7.84E8 




1115 


0, 


,13 


3.216 


3.56E9 


hd4.0818.1037 


989 


0, 


,04 


1.444 


6.68E9 


hd4.0298.0744 


782 


0, 


,54 


1.761 


4.21E9 


C4-06 


1357 


0, 


,06 


2.803 


I.IEIO 




1358 


1 


28 


2.803 


4.60E9 


hd4.1006.0680 


1063 


0, 


,07 


2.969 


4.16E8 


hd4.1486.0880 


1282 


0, 


,02 


-1.000 


0.00 


hd4.0367.0266 


813 


0, 


,04 


2.931 


7.34E9 




814 


1 


,06 


2.931 


1.8E10 


hd2.1928.1041 


109 


0, 


,23 


2.009 


2.81E9 


C2-05 


110 


0, 


,27 


2.005 


9.70E9 


hd2.1 739.1258 


162 


0, 


,68 


1.980 


1.37E9 




163 


0, 


41 


1.980 


2.3E10 


hd2.1398.1164 


274 


0, 


,25 


2.237 


4.75E9 


C3-02 


1179 


0, 


,23 


2.218 


1.25E9 


hd2.0698.1297 


553 


0, 


,35 


3.430 


2.09E9 


hd2.0705.1366 


516 


0, 


,13 


3.368 


1.77E9 


hd2.0529.1567 


616 


0, 


,08 


2.161 


2.50E8 


hd2.0624.1688 


561 


0, 


,10 


2.419 


8.61E8 


hd2.0725.1818 


503 


0, 


,16 


2.233 


1.47E9 


hd2.0743.1844 


503 


1 


,30 


2.233 


1.47E9 


hd3.1633.1909 


1499 


0, 


,48 


2.175 


4.88E8 




1506 


0, 


,99 


4.050 


1.36E9 


hd3.1824.1945 


1513 


0, 


,38 


2.050 


4.15E9 



^Targets arc listed in order of increasing RA, as in Table 

^Identificati on number(s) o f objec t(sj in HDF NICMOS Mo- 
saic catalog of |Papovich et al ] J2005I 1 with coordinates matching 
those of target galaxy. 

■^Separation in arcseconds between our target coordinates and 
those of matching HNM source. 

■^Redshift used by lPapovich et al | j2005| ) to fit stellar popula- 
tion model. A value of -1.000 indicates that no spectroscopic or 
photometric redshift was available. 

■^Stellar mass from fPapovich et al 1 120051) . A value of 0.00 in- 
dicates that no spectroscopic or photometric redshift was avail- 
able. 
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Fig. 1. — Color images of all 22 candidate and confirmed LBGs in our kinematic sample 
made from HST/WFPC2 -B450, Veoe, and /814 HDF images controlling blue, green, and red 
pixels, respectively. Targets are shown in order of increasing RA, as in Table [21 North is up, 
east is left, and each box is 6" on a side. 
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Fig. 2. — HST/WFPC2 Veoe images of all 22 candidate and confirmed LBGs in our kinematic 
sample. A schematic slitlet I'.'l x 6" long is shown in each case at the position angle of the 
observation, which is labeled in the lower left of each frame. Targets are shown in order of 
increasing RA, as in Table [2l North is up, east is left, and each box is 6" on a side. 
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Fig. 2. — Continued. 
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Fig. 3. — HST /814 images (left) and closeups of extended kinematic features in 2D spectra 
(right) of the 11 LBGs in seven systems showing such features. When a pair of targets is 
covered by a single slit, both targets are labeled in the image. Targets are shown in order 
of increasing RA, as in Table [2j Each image is 10" on a side and is oriented to show the slit 
vertical for comparison with the 2D spectrum. The slit width is I'.'l xcosz, where i is the 
relative position angle of the slitlet with respect to the slit mask. Each spectrum covers the 
same 10" vertically as the image, and roughly 250A horizontally, the exact value depending 
on the projected pixel scale and therefore, again, the position angle of the particular slitlet 
used relative to its parent slit mask. The beginning and ending wavelengths (A, lower left 
and right of spectrum) of each section of spectrum and the position angle of the slitlet (° 
east of north, upper left of spectrum) are labeled. The features include spatially resolved 
or separated pairs of absorption lines of Al II, Si II, and O I; and spatially or spectrally 
resolved or separated pairs of emission lines of Lya. The emission lines are usually the only 
relevant kinematic feature in their parent spectra; the absorption lines however are chosen 
as representative from among several. 
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Fig. 3. — Continued. 



-56- 




3500 4000 4500. 50QO SSQi! 6000 eSOO 7000 



40 
20 
O 

*P 



m 



60 
451 
SO 




C—l — 1 — i — 1— T— 1 — 1 — 1 — 1 L 1 1 f 11 1- 1 1 1 1 .1 1 J — 1 1 ! 1 

[ ^^^^^ 


\ -' ' > '1 ' ' ' ' 1 ; 


LJ 1 1 L_J 1 J 1 1 1 1 1 I 1 d 1 i ] d 1 1 1 1 1 1 1 1 1 1 1 1 L I 1 J 1 t \ i 



LJ 1 1 L_J 1 J 1 1 1 1 1 1 1 J_ 

3&00 4000 4500 


1 J 1 J 1 
5000 




5600 




6000 


i 1. 1 1 . 

6500 


1 I 1 I 
7000 


1 1 1 1 1 1 1 1 ' 1 1 I 1 [ : 








1 1 J 1 1 


' ' l-^ 


' 1 . i i. , 1 . I I 1. 1 . .1 >. i 




, , 1 , , 


. . f '-.. . 


.. . 1 


:> 

i.- . r;:^ 


1^ 


1 1 1 > 1 


B500 

1 1 1 1 < 


6000 

1 1 1 1 L 


1 1 'i~r-i- 






--L- -L-J J-. 



■-L L — L L - 



Mm sSsSbd i*s# %sii ?O0O 



( ^^^^^^ 




- 1 1 1 1 1 1 1 1 1 r 1 1 1 1 ,1 1 ,1 1 ,1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 t 1 1 1 


iiift m0 m& wm. 


m!^:0^. ' ^ ( ' ' ^ '1 ' ' ' ' 1 ' • * ■ i ' 


4. .1. 1 1 I I I P.'.j ^...j. 



3500 4000 4500 5000 5500 6000 6500 7O00 



Fig. 4. — Onc-dimcnsional spectra extracted from our two-dimensional LRIS spectra of all 
22 candidate and confirmed LBG kinematic targets. Spectra are shown in order of increasing 
redshift, except in the first panel, for which redshifts are unknown (redshift quality Qz = 1) 
and the objects are shown in order of increasing RA. Redshifts with Qz = 2 are followed 
by '?'. Each panel includes a spectrum of the night sky showing strong teluric emission 
lines. The wavelengths of common LBG emission and absorption lines are marked on each 
spectrum. The spectra have been cleaned of strong night sky emission line residuals and 
smoothed by 7 pixels (~ 8 A) to suppress noise. 
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Fig. 5. — Pairwise velocities vs. projected separation for close pairs of simulated LBGs wit h 
projected separation Ar < 7" for LBGs in Model 3 of Ivan Kampen fc Crawford I (120071 ). 
Open triangles indicate pairs with cosmological redshift difference A2; < 0.001, while dots 
indicate close pairs with larger redshift differences. The values for the four close pairs in our 
sample are shown with a x symbol. 
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Fig. 6. — Calculated dynamical mass (from Eq. [T]) vs. combine d halo mass for close pairs 
of LB Gs with projected separation Ar < 7" in the simulation of Ivan Kampen fc Crawford 
(120071 ). Symbols are as in Fig. [5l The solid hne represents Mnaio = ^Dyn- No strong 
correlation is seen. 
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Fig. 7. — Estimated dynamical mass vs. redshift for our 13 LBGs with kinematic signatures. 
Masses of close pair systems are shown as light blue filled squares at the pair's average red- 
shift, multiple knots as dark blue filled triangles, the galaxy (C4-06) with resolved absorption 
as a green filled circle, and masses derived from Lya linewidths as red filled pentagon s. Also 
plotted are the mean value for the stellar masses derived by iPapovich et al I (120011 ) using 
HST/WFPC2 + NICMOS optical an d NIR photom e try (o pen triangle); the median value 
for stellar masses of luminous LBGs by lShapley et al I (120011 ) using ground-based NIR (J and 
Kg) photometry [open square); the r nean value fo r the da rk matter halos of LBGs a.t z = 2.9 
derived using clustering analysis by lAdelberger I (l2005al ) {star symbol); the mean dynami- 
cal mass of BM/BX galaxies at 2; ~ 2 measured using redshifted Ha emission l ine widths 



and t he stellar mass using optical-MIR photometry including from Spitzer by lErb et al 



( 120061 ) (upper and lower crosses, respectively); the mea n dynamical mass of BM/BX galax- 
ies in the SINS survey measured with VLT/SINFONI by lForster Schreiber et al I (120061 ) {plus 
symbol); and the theor etical median and range predicted by the semi-analytical models of 



Somerville et al I (120011 ) {open circle). 
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Fig. 8. — (left) Stellar mass vs. estimated dynamical mass for individual LBGs with kine- 
matic signatures, (right) Combined stellar mass vs. estirnated s ystem dynamical mass for 
close-pair systems. Stellar masses are from iPapovich et al I (l2005l ). Symbols are as in Fig. U\ 
and the diagonal lines represent M* = Mdyn- Each point is labeled with the probability that 
the corresponding pair is a true physical association rather than a chance superposition, as 
estimated from the simulation results shown in Fig. [51 



